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Abstract 

We write the loop equations for the (3 two-matrix model, and we propose a topological 
recursion algorithm to solve them, order by order in a small parameter. We find that 
to leading order, the spectral curve is a "quantum" spectral curve, i.e. it is given by a 
differential operator (instead of an algebraic equation for the hermitian case). Here, 
we study the case where that quantum spectral curve is completely degenerate, it 
satisfies a Bethe ansatz, and the spectral curve is the Baxter TQ relation. 
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1 Introduction: Generalization of the two-matrix 
model 

Random matrix models have played a very important role in physics and mathematics 
[9l [26 | I36l438|. [50 | l62 | [67] . The eigenvalue statistics of large random matrices tend to 
universal statistical laws, which can be observed in many systems in nature, ranging 
from nuclear physics to finance and biology. Moreover, random matrices, treated as 
formal series, provide generating functions for counting discretized surfaces [T4"ll2"2"ll2"5| 
I2S1I3SG3], and are a very useful tool for string theory [30H33l[57H60l[6l] • 

So far, the random matrices studied were mostly hermitian matrices, because it 
was the easiest case. From the beginning, Wigner [62j|67J[69] introduced 3 ensembles 
of matrices, classified by an exponent 2/3, Hermitian: 2/3 = 2, real-symmetric: 2/3 = 1, 
and quaternionic: 2/3 = 4. However, it is rather easy to analytically extend the joint 
eigenvalue distribution to arbitrary values of /3 [T5 l l2"4" t l3"5] . 
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The case of one random matrix belonging to an arbitrary f3 ensemble, although less 
studied than the hermitian case, has received lot of attention, and many works were 
done. In particular loop equations have been known for a long time, and their recursive 
solution was recently proposed in pT7HT9| [^6] . 

For the hermitian case, it turned out that a 1-matrix model was not general enough, 
for instance the spectral curve for 1-matrix model is always quadratic, it can't have 
higher degree, the critical behaviors are related to only a small subset of all possible 
conformal field theories. A 2-matrix hermitian model turns out to be much more 
general, it allows to reach spectral curves of any degrees, and any conformal minimal 
model 56J. Moreover, as applications to counting discretized surfaces, a 

2-matrix model allows to count discrete surfaces carrying colors [25|l54"j. 

Unfortunately, trying to generalize the 2-matrix model to arbitrary /3 seemed very 
difficult, and was almost never studied. One reason, is that the joint eigenvalue dis- 
tribution was not known, because the integration over the angular parts of the matrix 
was not known (properties of angular integrals for hermitian matrix models can be 
found for example in [68]). 

Recently, some progress was made [6] in those angular integrals, and although we 
don't know yet how to compute angular integrals for arbitrary /3, we know already 
that those angular integrals have to satisfy some differential equations, and this is 
sufficient to derive loop equations. This is what we do in this article, we derive the 
loop equations, and solve them perturbatively by expanding in some small parameter. 



1.1 The hermitian two-matrix model 

The hermitian two-matrix model is given by the partition function: 

Z Herm = e F = J d M 1 dM 2 e-^ tI ^ M ^ +v ^- M ^ (1.1) 

where Mi and M 2 are N x N hermitian matrices and the measure dMi is the corre- 
sponding Lebesgue measure associated with all independent entries of the matrices. 
V\ and V 2 are called the potentials, we shall assume in this article that V\ and V 2 ar e 
polynomials. The parameter T, often called the temperature, is redundant, it can be 
absorbed by a change of variable M 2 — > TM 2 and a redefinition of V 2 , however, we pre- 
fer to keep it for convenience, and because of its future geometric interpretation. The 
hermitian 2-matrix model was introduced in particular as a formal series to study the 
Ising model on a random surface, i.e. the Ising model coupled to 2D gravity [26 l l4*4" t l54"] . 

It is very well known [61] now that this integral can be rewritten in term of its 
eigenvalues problem as: 



Z Herm = / dXdYA(X) 2 A(Y) 2 e 



N N 

£fi(*i)+E^(w) 



I(X,Y) (1.2) 
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where X = diag(xi, . . . ,xn) and Y = diag(yi, . . . , i/n) are diagonal matrices rep- 
resenting the eigenvalues of Mi and M2, A(X) = Yli 

the Vandermonde 

i<j 

determinant, and I(X, Y) is the Itzykson-Zuber integral [5T||52] which corresponds to 
the integration about the angular variables: 

I(X,Y)= [ dUe^^ XUYU ^ (1.3) 
Ju N 

where is the unitary group equipped with the Haar measure. Relatively to the 
Itzykson-Zuber integral, one can also define the quantities: 



[ dU\\U hJ \\ 2 e^ T ^ XUYU ^ (1.4) 
Ju N 

which can be used to determine the Itzykson-Zuber integral by the formula: 

N N 

i(x,y) = j2 m ^=J2 m ^ ( L5 ) 
1=1 3=1 

Note in particular that the r.h.s. does not depend on the second index i or j. Tech- 

N N 

nically, the last formulas are obvious since ^||K,j|| 2 = 1 = Sll^j'll 2 m the unitary 

i=i ' 3=1 

group. Eventually, in the hermitian case, it is also known that the Mjj's satisfy the 
Dunkl equation [IHlElHlEI] , namely: 

VI < ij < N : Am, + E ^ I f * ' = ViMv (1.6) 

UX{ X{ Xfc 

We will see in the next sections that these properties can be extended uniquely in 
the case of the /3- deformation of the hermitian two-matrix model case. 

1.2 Generalization to arbitrary (5 two-matrix models 

Similarly to what happens in the one-matrix model, we would like to generalize the 
hermitian two-matrix model to other /3-ensembles of matrices. As in the case of one- 
matrix models, the usual generalization is made from the diagonalized version of the 
problem ( 11.21) by: 

„„ r jv n 

z/= dXdYA{X) 2 ^A{Y) 2 ^e ^ ^p{X,Y) (1.7) 

where Ip(X,Y) is a generalized version of the Itzykson-Zuber integral that we will 
describe later. Note that with our convention, the hermitian case corresponds to (3 = 1 
(notation used by people working on AGT conjecture [2], and Laughlin wave function), 
whereas sometimes in the literature it is normalized to /3 = 2 (Wigner's notation 
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j). We prefer to use the first notation in order to avoid unnecessary powers of 2 
in all formulas. 

The diagonalized version (ll.7p corresponds to the idea of the generalized matrix 
integral: 

V' = " / rfM lC iM 2 e-^ tr[v/l(Ml)+v/2(M2) - MlM2] (1.8) 

where E^,p would be an ensemble of matrices such that the diagonalization of Mi and 
M 2 gives fll.7p . Note that this definition is only formal since apart from /3 = 1, 0, | 
or 2 where En corresponds to hermitian, diagonal, real-symmetric or quaternionic 
matrices, no other ensemble is known presently (it would be interesting to see if the 
/3-matrix ensemble introduced by [35] also reproduces the Ip term). 

1.3 The angular integral 

We shall now define the Ip{X, Y) angular integral, so that it coincides with the actual 
angular integral for the matrix cases (3 = 1, 0, |, 2, and in fact we shall use the angular 
matrix integral defined in [6]. 

Following [6], we first define a good generalization of the M it j and from (I1.5P we 
will define the generalized version of the Itzykson-Zuber integral. In [6], the authors 
claim that the natural generalization of the Mij, noted here M^j is the following: 

1. The MjjS satisfy the generalized Dunkl-Calogero-Moser equations: 

V 1 < i,j < N : A M g) + £ = "frtf (1.9) 

dxi ' J Xi - x k T J 



N/3 
T 

AT/?! 
T 



(The factor 4^ comes from the same factor in the exponential term of the parti- 



tion function which can be absorbed by the change Y -H- -pjrY) 

N 

E 

i=l 

N 

and Ip = E ^fj must be independent from i 



2. Mjfj must be stochastic matrices, i.e. Ip = J2M^ must be independent from j 



3. They must have the symmetry M^(X,Y) = M^-j (Y, X). In particular, this 

N [0), 

implies that Ip(X, Y) = ^2M^ (X, Y) must be symmetric in the exchange X -h- 

i=i 

Y. 

4. They must also have the symmetry ^(X a ,Y) = (X : Y) where a £ Tin 
is a permutation. This implies that Ip(X a , Y) = Ip(X, Y). 



5. They are normalized so that Ip(ld.NxN, Y) = e 



^ try 
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These conditions define the M^ s uniquely making the previous set of conditions a 
proper definition (in [5J, an explicit solution for those is provided as an integral, 
or is given by a recursion on N). Moreover all these properties are standard results in 
the hermitian case P = 1 and hold for /3 = 1, §, 2. In particular note that M$(X, Y) = 
Mj^(Y, X) gives another formulation of (jl.9p as: 

V 1 < i,j < N : A M ff + p-£ M V- M « = ^ X Mf (1-10) 

Similarly to what happens in the hermitian case, it is logical to define the generalized 
Itzykson-Zuber integral by: 

N N 

ip(x,Y) d = f (x,y) = J2 M S ] ( x >n ( L n) 
i=i i=i 

Note again that this definition recovers the known cases when /3 = 1, |, 2. Eventually, 
from this definition we can prove (see [6j) that the generalized Itzykson-Zuber integral 
Ip(X,Y) satisfies the following Calogero-Moser equation [IE] : 

where -ffjf ^ s Calogero-Moser Hamiltonian. Note again that Ip(X, Y) is well known 
to satisfy f)1.12p in the (3 = 1, |, 2 cases. 

Now that we have introduced a proper generalized two-matrix model, we can try 
to solve it in the large N limit and in various regimes of the parameters. In particular, 
since no (bi)-orthogonal polynomials techniques are known in the case of an arbitrary 
exponent /3 in the Vandermonde determinant (although some work has been done 
in [31]), we will use in this article the method of the loop equations to determine the 
correlation functions of the models. 



2 Summary of the main results 

Our goal is to compute the asymptotic (possibly formal) expansion of the 2-matrix 
model integral and its correlation functions. 

V = " / dM 1 dM 2 e" ^ tr m^+v^-MM] 
and the cumulants (subscript c ) of expectation values of products of traces of resolvents 

W n (xi, . . . , x n ) — ( Tr — — . . . — — ) . 

V ^ \ x 1 -M 1 x n -M 1 / c 
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We want to expand them in powers of a small parameter e as 

9 

9 

The small parameter e is a combination of and /3. There can be several asymptotic 
regimes of interest, which we discuss in detail in section HI 

Main results 

• In section |3] we derive the loop equations for the 2-matrix model with arbitrary 

P. 

• In section |5] we solve the loop equation to leading order for the resolvent Wi, i.e. 
we compute Wi°\ We show that 

r(0)/_ ^ T 



W?\x) = - 



N ip{x) 

where ip{x) is solution of a linear differential equation with polynomial coefficients 

where £ (x, y) is a polynomial in 2 variables. 

• This equation for ip can be thought of as a " quantum spectral curve" . We remind 
that for — 1 one gets an algebraic equation S(x, W±°\x)) = 0. Here, we replace this 
algebraic equation by a linear operator S(x,y) where 

, T d r „ . T 

i.e. the non- commutative version of the algebraic equation S(x,y) = 0. 

• This equation can also be seen as a Baxter T-Q equation where ip(x) is the 
Q function. Its zeros Sj's (ip(si) = 0) are thus the Bethe roots, and the equation 
S(x,y).tp = can be translated into a Bethe Ansatz equation for the s^'s. We do it 
in 16.171 The Bethe equation for S = diag(si, . . . , sn) can be written in the following 
way: find a matrix B such that: 

M(B)-S)e = , e t (Vl(S)-B) = , [S, B] = 1 (Id - e e*) 

where e = (1, 1, ... , 1)*. 

• In section [6761 we show that this Bethe ansatz equation can also be written as the 
extremization of a Yang- Yang action: 

^S(S,S,A,u) = tr Vi(S) + tr V 2 (S) - tr (SASA -1 ) - ln(A(5')) 



S 



~ In(A(5)) + I In det(A) - ^ (Ae - e) 

At the extremum of 5 we have S = diag(si, . . . , sjv). 

The leading order Fo of In Z is the value of the Yang- Yang action at its extremum, 
and the first subleading term i*\ is half the log of the determinant of the Hessian of S: 

In Z w = ^ e 29_2 F s , F = -S(extremum) , F x = - - In det 5" 

• In section 16.7.11 we compute the leading order of the 2-point function W 2 in 
terms of the Hessian 7i = S" of the Yang- Yang action: 



N 



2 („l _ .\2- 



x — Si r [x — s 



In a similar way, we compute the leading order of the 3-point function in terms 
of H = S" and S'" in f l6T^|) . 

• In section [7] we provide a "topological recursion" algorithm to compute every 
subleading correction to any ra-point function, i.e. Wn \ in the form: 



N 



n 



W" y > ( Xq, X\, . . . , x n ) — N Res K{xq, x) ( combination of Wffl s with 2h-\-m < 2g 

i=l 

However, in contrast with [IB], the recursion kernel K(xq,x) is here a matrix of di- 
mensions g?2 = deg V^. We derive the topological recursion in section 17.41 and then we 
describe the kernel K(xq,x) in more details in section 1731 and in appendix. 



3 Loop equations for the /3- deformed two-matrix 
model 

The loop equations, also called "Schwinger-Dyson" equations, are a very powerful tool 
to study random matrices, perturbatively in the expansion in some small parameter 
(for the hermitian case, the small parameter is usually 1/N in the large N limit), 
see [3j[23j|2|2llQj[53j|66]- In short, Schwinger-Dyson equations exploit the fact that an 
integral is invariant by change of variable, or an alternative way to obtain the same 
equations is doing integration by parts. 



3.1 Notations for correlation functions 

For the hermitian two-matrix models, the loop equations are well known p0ll4Tl - l4"3ll4~7l 
HEIISSIES] and they are very useful to compute the large N expansion of correlations 
functions, wlf 1 and symplectic invariants F g presented in a more general context in [4"8] . 
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In the one-matrix model, it is also known [36] that the loop equations method can be 
generalized for arbitrary (3 quite directly from the hermitian case. Indeed, we remind 
the reader that we can see loop equations as integration by part or Schwinger- Dyson 
equations of the diagonalized problem (11.71) . In the present /3-deformed two-matrix 
model case, we will see that we can get also some loop equations by following the same 
approach. Before writing down these equations, we need to introduce some notations 
(the same as in [20|l4Tll42]): 

• The potentials are assumed to be polynomials: 

V((x)=J2 t * xk > V^x) = Y,ikX k (3.1) 

k=0 k=0 

• The correlation functions are defined by: 

W n (z 1 ,...,z n ) = ( —^-—Z ) ( 3 ' 2 ) 

\w,...,in=l Zl Xh Zn Xin / C 

where the brackets ( ) indicates that we take the expectation value relatively to 
the measure defined by (11.71) for the random variables X = diag(xi, . . . , x n ) and 
Y = diag(yi, . . . , y N , namely: 



(A(X,Y)) d ^±- [ dXdYA(X,Y)e-^ tTV ^ x)+v ^ Y)) A(X)^A(Y) 2l3 I p (X,Y) 
Zh J 



(3.3) 

The index c stands for the connected component (also called cumulant) of the 
correlation function that is to say: 

< A x > = < A x > c 
< A X A 2 > = < A X A 2 > c + < Ax > c < A 2 > c 
< AxA 2 A 3 > = < AxA 2 A 3 > c + < A X A 2 > c < A 3 > + < A X A 3 > c < A 2 > 

+ < A 2 A 3 > c < A 3 > + < Ax >< A 2 >< A 3 > 

n k 

<Ax...A n > = < Aj >= II < Al > >c 

k=X IitiU 2 ---Vlk=J i=l 

(3.4) 

In order to have more compact notations, we will sometimes simply denote 
W(x) = Wx{x) for the first correlation function. 

We can define similarly the second type of correlation functions: 

E ■■■ — - — ) ( 3 - 5 ) 



and denote W(y) = Wi(y) for the first correlation function. 
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To close the loop equations we will need to introduce the following functions: 

U n (x,y;zi,...,z n ) = 

y I 1 M%>v>(y)-v>( yj ) 1 
m i-i\ x ~ Xi y ~ y i z l -x il "' z n -x in 

(3.6) 

Note that they are polynomials in the variable y, and remember that M^(X, Y) 
depends on X and Y and is also a random variable. In a similar way, we also 
introduce: 

£ / V{{x) - V>{ Xi ) M%> VM - V±{y 3 



(3.7) 

Note that this time P n (x, y;z±, . . . , z n ) are polynomials in both x and y. 

The loop equations method is a powerful method of deriving an infinite set of 
equations connecting the correlation functions and the functions P n and U n introduced 
before. There are several ways to derive the loop equations and we choose here to use 
the integration by part way which is the most convenient in our setting. In the case 
of hermitian, real-symmetric or quaternionic matrix models, these loop equations are 
already known and can be derived directly from the matrix integral settings (before 
any diagonalization). They can be found in many different places in the literature 
[T^Mini[MgTIB3lH8ll5lll56l[65] . In our arbitrary (3 setting, we do not have the initial 
matrix integral model properly defined, thus we must adapt the derivation of the loop 
equations with our definitions (jl.9p and (11. lip . This is done in different successive 
steps that are fully detailed in the next sections. 

3.2 First step: deriving an auxiliary result 

We start by evaluating the following integral: 

= Y [ dXdY-^- ( e-^^^W+^AiX^AiY^—^—M^) (3.8) 

^ ./ OVi V X - Xi ,J J 



The result is clearly null if we choose an integration path which goes to oo in directions 
in which e - "^ 1 ^ and e~^ V2( - Y "> vanish exponentially, because we integrate a total 
derivative. The right hand side yields 3 different contributions: 

• Acting on the exponential we find: 



T — , \ X — Xi 1/3 

1,3=1 \ 
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Acting on the Vandermonde determinant, we find: 



1 M 



i.j 



— 1 % -VkX- X % If, 



(3.10) 



Eventually, acting on and using the Dunkl differential equation satisfied by 



the M/f's QUOD we find: 



N 
\i,j=l 



x - Xi \ T ^ Vj -y k J I B 



(3-11) 



Then we see that (13.101) cancels with the last part of (13. lip so that we have our 
first equation: 



N 



E ( n^ M 4) -t I * MS> 



'■7 = 1 



' 1 1,3=1 x 



Remembering now that = Ip gives: 




(3.12) 



(3.13) 



3.3 Second step: finding the loop equations 

Then consider: 



= Y [ dXdY-^-( e-fM^WWlAm^m^I 
4-1 J dxi\ x - Xi 1 



whose r.h.s. yields 4 contributions: 

• Acting on the exponential, we find: 



N(3 £ lv[{ Xi ) Mff V>(y) - V>( yj ) 



T ^ \X-Xi Irs 
i,j=l \ 



y-Vi 



(0 



V-Vi 



(3.14) 



(3.15) 



Acting on the Vandermonde determinant we find: 



^e/e- 2 



1 M ™ - Vi(y,) 



X{ Xfc X X{ 1 a 



11 1 



V-Vi 



(3.16) 
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Acting on — we find: 



Acting on and using (II. 9p we find: 



(3.18) 



Now we observe the following identities: first in ( 13 . 1 5[) we can perform V{(xj) = V{(xi) — 
V{(x) + V{(x) so that we have (using the notations Uq and Pq introduced in (13. 6p and 

my- 

(i)e>-^f (V;(x)U (x, y) - P (x, y)) (3.19) 
Secondly we can split (ii) in the following form: 

1 Mm VM - VJfe) 



(«) ^E(E^ 



N I i i uWv/ 



i ^'(y) - ^'(%) 



1 Mff n(y) - %) 



iJ=1 •'•-•'■/■• •'•--•'•/.. /i // .'/, 



E^ 



+ ^ £ / E Wte)-W > (mT 



(3.20) 



i fa - x k )(x - Xi) y-yj 



Note that (ii)" is the same as the last terms of (iv) so that it cancels out. We can also 
split (ii)' into a sum over i, k minus the case i = k which is identical to (Hi) except as 
a factor 3. Therefore we can regroup (ii), (Hi) and (iv) to get: 

(a) + («o + ,u„ = (i - « t li^mzmA) (1) 

^ \ ( x - x i) h v- Vi I 
i~^=i\ x - x i x - x k b v-Vi J 

W £ /^MlYM^YM) ,3) (3.21) 
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Observe now, that we have: 

Then we have also: 

T I 9 . A = f /J 1 ^ ^(if)-%) 

_ y / 1 < } V^(y) - ^(yj) \ / 1 < } 

.^^i-SiX-ii I/j y-y^ 
-wi>)[/ (x,y) (3.23) 

We recognize here one the second term of (13.211) : 

(2) # ^(g^Uo(x,y)y + (3W(x)U (x,y) 

& 8U 1 (x,y;x) + 8W(x)U (x,y) (3.24) 

Eventually, we are left with (3). We perform the change yj -H- yj — y + y and split it 
into two parts: 

3 ^ -JLyUfay) - -fViMY, { T~ 

T T ~^ \ X ~ X i b 

But remember that ^2M-\j = I p. and that from our first step we have (13. 1 3f) so 
that eventually: 

(3) O ^W*, 2/) - + ^-(-N + xW{x)) (3.26) 

Now we can put everything back together to get the loop equation: 

NB d 
= ~ (V((x)U (x, y) - P Q (x, y))-(l- 0)g^U o (x, y) 

NB X 

+BUx(x, y; x) + BW(x)U (x, y) + —yU (x, y) 

NB „ , , , NB, 
~Vi{y)W{x) + ^f(-N + xW(x)) 

(3.27) 
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which can be rewritten (performing a multiplication by — jfp) as the master loop 
equation: 



y - V((x) + ^W(x) + hd x ) U (x, y) = (V^(y) - x)W(x) - P {x, y) 
+ N-^U 1 (x,y;x) 



(3.28) 



where we have defined h = -^(1 — ^) which we prefer to write 




(3.29) 



3.4 Higher order loop equations 

With the same approach as before, one can deduce higher order loop equations for 
U n (x, y, £1, . . . , , £ n ). Indeed, if one look at the auxiliary integral: 



N 
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£ l dXdY 7^[ e 
1 



— (p-^ 



1 



A(X) 2 ^A(Y) 2 ^—^—M-j(X, Y) 



n 



. tr 



C P (^,...,QA(X)^A(Y)^- 



1 



-M 



1 -j 



(3.30) 



with 



C P (6,..-,e P ) = / dXdZe~^^ v ^ v ^A(XrA(Zri,(X,Z) (f[ tr ^J-_ j 



One finds: 

A' 



1 m; 



f (/3) t> , 



xW n+ i(x,fi, . . . ,£„ 



(3.31) 



(3.32) 



Then the second step can be carried out using the integral: 



N r 8 

4,7 = 1 



tr (Vi(x)+y 2 (y)) A(X) 2/3 A(y) 2/3 



(3.33) 
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n 

First, we find the same terms as before (the new term Y[ t r g _ x being treated as a 

P =i p 

constant). This gives: 



T 



(v((x)U n (x, y,£i,..., - P n (x, y, £1, . . . , 



d 

-(i - 0)-^u n (x, £i> • • • ' £») + /^n+ifo 2/; a? > • • • > x O 

+/3W(x)C/ n (x, y, £i, . . . , f n ) + —yU n (x, y,^,..., f n ) 
+f3 Yl Wi + \i\(x, I) U\j\(x, y; J) 

IUJ={^,...,^},J^0 



T 

(3.34) 



■^'(j/jWWi^, 6, • • • , £n) + sWn+l^, £l> • • • , fn)) 



n 

However, the main difference is that now the derivative will also act on TT tr g 1 „ . 

P =i 



This gives an additional contribution: 



i i 



Then observe the two identities: 

1 d 1 



and: 

11 1/1 1 



X X{ Xi X \s& X Xi 

Therefore the action of the derivative can be rewritten as: 



(3.35) 



(3.36) 



(3.37) 



V-^ 9 ( f/ n _i(x, J/, ^i, ... , Ck-U ■ ■ ■ i Cn) ~~ ^n-l(6fc> 2/) £l> ■ ■ ■ ; £fc-l; £jfc+l> • • • > £n) 



fc=l 



(3.38) 

Eventually, adding (13.341) and (I3.38f) to get (I3.33f) we get the higher order loop 
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equations: 



T \ T 

V ~ V(( x ) + jj w ( x ) + hdx J U n (x, y; . . . , f n ) + — U n+1 (x, y; x, £ n 

| Y, W 1+lIl (xJ)U lJl (x,y ] J) + ^-J2-^- 
/uj={a,.-,«n},/^0 k=i Sfe 

[ x , y, £i> ■ • • , Ck-ij Cfc+i, • • • j £,n) — U n _i((,k, y, £i> ■ ■ ■ > ■Cfc-i, C^+i, ■ ■ ■ , Cn 

z - £ fc 

X i Vi Sl) ■ ■ ■ i 4"n) 



(3.39) 

These loop equations (I3.28p . (13.391) are exact equations satisfied by the correlation 
functions, they hold for every N and every (3. Now, in order to solve them, we want 
to consider some limit and expansion around that limit, in which they can be solved 
order by order. 

3.5 Relationships between correlators 

The loop equations which we have written, involve the resolvent W (x), as well as higher 
correlation functions W n , and also some auxiliary functions U n and P n . However all of 
them are related. 

First, notice that P n is the large x polynomial part of V{U n : 

P n (x, y;xi,..., x n ) = (V((x) U n (x, y;xi,..., x n )) + (3.40) 

Second, notice that the large y behavior of U n is related to W n+ i. 

U n (x, y; Xl ,..., x n ) ~ t d2 y**- 1 W n+1 (x, x u . . . , x n ) + 0(^ 2 " 2 ) (3.41) 

y— >oo 

And then, the loop insertion operators allow to increase the index n by 1. Define: 

d ^1,0 



dVAx) f-f x k+1 dt k - X 



y- 



k- 



k = l 



d A 1 . d 



SV 2 (x) f^y"* 1 dt„^ 

These operators, called loop-insertion-operators, are formal and they have the property 
that: 

9_YM _ s 1 9 _YM= s 1 ^ 

dVi{x') hl x-x> ' dViix')' hl {x-x') 2 K ' 
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The main interest of these operators is that they can be used to connect the correlation 
functions W n (xi, . . . , x n ) to the next correlation functions W n +i(xi, ■ ■ ■ , %n+i)' 

W\x) 



T y ' dVtix) 

— 1 \X\ • • • j ^n— lj 

W n {x x , ...,x n ) 



We also have 



T dVl{Xn) 

N/3 TT t . dU n _i(x,y; x 1 ...,x n - 1 ) 

U n {x, y,xi,..., x n . 
T oVi(x n ) 

(3.44) 



N/3 , N dP n „ 1 {x,y;xi...,x n -i 

*n\Xi U) X\i . . . , X< 



— i\Xj Hi X\ . . . , x n —\ —i\X n , Hi X\ . . . , x n _x) 



(3.45) 



4 Topological expansion (s) 

One sees, that the loop equation (I3.28P involves U and U\, and similarly (I3.39P involves 
U n and U n+ i. This is why loop equations can only be solved perturbatively, order by 
order in some small parameter, in a regime where the U n+ i term is subleading compared 
to the U n term. 

We thus need to expand all our observables in powers of some parameter, and this 
expansion is usually called topological expansion. Most often for hermitian matrices, 
the small parameter is chosen to be l/N where iV is the size of the matrices, i.e. the 
topological expansion is an asymptotic expansion for large random matrices. 

Such a power series expansion does not always exist for convergent matrix models, 
but by definition it does for formal matrix models which is our context here (we do not 
interest ourselves in convergence here). For the subtle differences between convergent 
and formal matrix models in the hermitian case, we invite the interested reader to refer 

to HUES]. 

We thus need a small parameter in which to compute an expansion. For the her- 
mitian case, the good choice is an expansion into powers of l/N in the large N limit. 
However, for ^ 1, we can already see on loop equations, that the most interest- 
ing regime (in which loop equation are most conveniently solved) is a regime where 
h = 0(1), and so is not simply N — > oo. In fact, depending on which application of 
random matrices one is interested in, several limit regimes can be interesting. Fortu- 
nately, all of them are related together, and knowing one expansion allows to recover 
the others. Let us present 3 main regimes: 
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Before presenting those regimes, let us introduce some usual notations: 

g s = string coupling constant 



9s - 


1 




ei = 






9s 


£2 = 


~ VP 


h = 


Cl + ^2 


9 2 s = 


= -ei£2 










£2 



1 st equivariant parameter 
2 nd equivariant parameter 



(4.1) 

Those notations are those used in different applications of matrix models, in particular 
string theory piIJ[2Hl[2H] and in relationship to the now famous AGT conjecture [2J. 

4.1 Topological regime 

This is the regime where we expand in powers of g s at fixed h. This regime is the 
one interesting for applications to geometry and topological strings. In this regime, we 
wish to expand: 

oo oo oo 

F = In Z = 9 2 r 2 F 9 (h) = J2J2 Si 9 " 2 hk (4-2) 

9=0 g=0 k=0 

Each coefficient F g (h) = F gk have a geometric interpretation. Notice that F g ^ 

depends neither on N nor on (3 (but of course it depends on the potentials V±, V2 and 
on T). 

In topological strings, F g j. is the generating function for counting equivariant 
Gromov-Witten numbers IT21. 



4.2 Large N regime 

This is an expansion in powers of N, at fixed /3. This is the regime interesting for large 
random matrices. Most often, one is interested in specific values of f3, like /3 = 1/2 or 
(3 = 2. We shall write: 

00 

F = InZ = J2( N / T ) 2 ~ l F i{P) ( 4 - 3 ) 
;=o 

This expansion is related to the topological expansion by 

E(— ) 

FiW) = E P 1+9 ~ l (P - l)^ 2s ^-2 9 - (4.4) 

9 =0 

For example 

F ((3) = (3 F 0t0 , F 1 (J3) = (J3-l)F 0tl (4.5) 
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(/3-D 2 - 0-1 (/3-D 3 

P P P 2 

and more generally by (I4.4p 

4.3 WKB large j3 regime 

This is an expansion in powers of at fixed N. This is the regime we shall consider 
below. It can be computed with standard WKB techniques, but we shall analyze it 
here with loop equations. 

We expand in powers of g s = at fixed iV: 

oo 

F = lnZ = 5>f- 2 A(iV) (4.7) 

h=0 

This regime can be related to the topological regime as follows: 
For example 



(T \ k TNT 
hk = F (T/N) , f 1 (N) = F 1 (-)- ¥ F>(- 



4.4 Notations for the expansion 

In this article, we want to study perturbatively the BKW regime where — > oo, N 

T 



fixed, and we write g s = jfrs- We thus expand all our correlation functions as: 



W n (x u ...,x n ) = [-) ^g: r^ n W^{x u 

9=0 



Uo(x, y) - (u^(x, y)-x + V^y)) + - £ g? U^(x, y) 



9=1 

/Arxn+1 oo 



9=0 
n+1 oo 



P n (x, y; an, . . . , x n ) ( - ) ^ P ^ x ' 2/; • • • , O 

^ ' 9=0 

(4.9) 

Note that the functions Wn , Un \ and Pn are only functions of iV (and 
of course of their corresponding X;, of T and of the potentials V±, V 2 ) but 
not of g s , or H. 
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Eventually, we define the numbers f g by the expansion of the logarithm of the 
partition function itself: 

oo 

Z = e F , F = Y,9 2 r 2 f 9 (N). (4.10) 

9=0 

Again, f g (N) is also implicitly a function of T, V\ and V2, and we recall that f g {N) is 
related to F g (h) which is the usual topological expansion. 

As we will see in the next sections, most of the formulas get nicer if we shift the first 
correlation function by the derivative of the potentials. Therefore, for convenience, we 
introduce the shifted functions: 

Y(x) := (Vl(x) - W[°\x)) , X(y) := {V}(y) - W?\y)) (4.11) 

It is also useful to redefine a shifted version of the polynomial P(°\x,y) as: 

E(x, y) := (V;(x) - y)(V^y) - x) - P (0) (a:, y) + T + ^, (4.12) 
which is a polynomial in both x and y (and it depends on T, N, Vi, Vi but not on pi). 



5 Spectral curve and Loop equations as an ODE 

In this section, we shall see that the loop equation satisfied by the leading order 
of the resolvent Wi(x) = ^tr x _} Mi J can be written as an ODE 

where S(x,y) is up to some shifts, the polynomial E(x,y) introduced above in f)4.12p . 
This ODE is called the "spectral curve" and plays a central role. 

5.1 Loop equations in g s expansion 

The leading order at large /3 and fixed of loop equation (I3.28j) (notice that we have 
W{ (x) = V{(x) — Y(x) and Uq (x,y) includes a shift by x — V^iyj) is: 



V - Y(x) + ^d^ U^ix, y) = (V;(x) - y)(V^y) - x) - P (0) (x, y) + T + ^ (5.2) 



which is equivalent to (See (]4.12p ): 




(5.3) 
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This is the master loop equation, also called the "spectral curve". It generalizes the 
hermitian 2-matrix models's spectral curve [X3 | l¥T | [5 ^[56l[65] . 



If we expand in powers of y: 

d 2 



d 2 +l 



X) 



k=0 



k=0 



equation (15. 2p reads: 



C<SLi(*) + ( -Y{x) + -d x ) U 0<k 



T 



(°), 



X 



EJx) 



Remark : equation ( 15. 5 P with k = d 2 + 1 gives us that: 



~td 2 



In the same way, since PQ°\x,y) is only a polynomial in y of degree d 2 — 1 we have 
also: 



(5.4) 



(5.5) 



(5.6) 



0gU(z) = i d2 (Vl(x) - Y(x)) - t d2 ^ 



(5.7) 



5.2 The linear differential system 

In order find the lower coefficients Uq°1(x) with = 0, . . . , ^ - 2 we introduce for 
convenience the functions ipk( x ) (0 < k < ofo) defined by: 



TT ( ,&cii)k{x) (o), » def Tip'(x) T d . , 



ip(x) 

which is equivalent to say that: 

N rz w (0) 



N i>(x) N dx 



(5.8) 



(5.9) 



By definition, the function %p(x) is only determined up to a global multiplicative con- 
stant which can be fixed by specifying the lower bound of the integral in the last 
formula. 



The first loop equation ( 15. 51) can now be rewritten into a linear differential equation 
(ODE): 



T 



4 )k _ 1 (x)+(^d x -V{(x))^\x) 



E k (x)^(x) 



(5.10) 
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and using i^>d 2 (x) = —td 2 ip(x), this leads to a set of linear equations that can be written 
into a matrix form: (equivalent to f)5.10p ) 



/ 



LA. 

Ndx 



\ ^o(z) / 



l d 2 -l 



J d 2 -\ 



(X) 



Ei(x) 
Eg(x) 



-1 

V{{x) 










C(x) 



d 2 

( i) d2 {x) \ 

^d 2 -i(x) 



o 



-1 

V{(x) 



\ 



( i^d 2 (x) \ 

Tpch-x(x) 

\ Mx) J 



(5.11) 



Therefore, we have obtained a system of first order linear ordinary differential 
equations given by a companion-like matrix. The determinant of the matrix 
C(x) is given by the standard formula for the determinant of a companion matrix: 



-t d2 det((y - V((x))Id + C(x)) = E(x, y) 



(5.12) 



That is to say that the characteristic polynomial of the matrix C(x) is exactly our 
function E(x, y). 

There is no known notion of integrable system associated to the arbitrary two 
matrix model, but the fact that the spectral curve turns out to be a characteristic 
polynomial (also called spectral determinant), is very suggestive from the point of 
view of integrability jUtTJIH], this could be a promising route, and needs to be further 
investigated... 



5.3 ODE satisfied by ifj(x) 



The matrix equation 05. 11 j) can be used to find a linear ODE satisfied by the function 
ip(x). Indeed we can rewrite it in following way: 



V((x) -j^d x j ip d2 -i(x) 



v;( x )-L dx ) Mx) 

V{(x) - ^d x ) Mx) 



— = WA X ) + Wd2-i{x) 

td 2 

E d2 -i{x) 



td 2 

Ex{x) 
Epjx) 

td-2 



-lpd-l(x) + 1pd 2 -2(x) 



ipd 2 (x) +ipo(x) 

Tpd 2 (x) 
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(5.13) 

Multiplying line k on the left by (V[(x) — hd x ) d2 ~ k and summing the whole set of 
equations give that: 

(rp \ d,2 + l d 2 
V{{x) - -d x iftz) - (V!(x) - hd x ) k (E k (x)iP(x)) = (5.14) 
' k=0 

which can be rewritten as: 




(5.15) 

Note in particular that since x and y do not commute ([y, x] = S), we must specify 
the position of the y compared to x. In the last formula, it is implicitly assumed that 
powers of V[ — y are always to the left of all powers of x, and: 

(Vi(x) - yf = (Vl(x) - y) (V((x) - y) . . . {V{{x) - y) (5.16) 

V v ' 

k times 

Eventually we find that ip(x) must satisfy a linear ODE of order di + 1, with 
polynomial coefficients, given by the quantized spectral curve E(x, y). This generalizes 
the Schrodinger equation in the case of the arbitrary-/? one- matrix model |46j . 

5.4 Quantum spectral curve 

The spectral curve E(x, Y(x)) = in the hermitian case, gets replaced by 

%A = jj » E(x,v;(x)-y)-^(x) = o- (5.17) 

This can be interpreted as having now a quantum spectral curve. 



5.5 Non-unicity of the solutions of the loop equations 

In order to identify a solution of loop equations, it remains to compute the polynomial 
coefficients Ek(x) of this ODE, and choose a solution ip. 

Then, as soon as ip(x) is chosen, we can find ipd 2 ( x ) by using ^d 2 (x) = —td 2 il){x). 
The other functions ipk{%) follow by a descending recursion with the help of f)5.10p . 
Eventually, the connection between ipk( x ) an d uj?\x) is given by (15. 8p so that the 
polynomial U^°'(x,y) is known. 

One may think that the knowledge of the loop equations would be enough to deter- 
mine all the correlation functions by solving them properly. But we remind the reader 
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here that in general the loop equations (even in the hermitian case) have infinitely 
many solutions. 

Even if (15. lip gives us some linear differential equations with polynomial coef- 
ficients, the main problem, also present in the hermitian case, is that the function 
E(x,y) is not known at the moment. Indeed, in the definition of E(x,y), we face the 
polynomial P^ (x, y) whose coefficients are unknown. The coefficients of the polyno- 
mial PQ°\x,y) are not given by the loop equations, they have to be determined by 
some other informations which we have not used so far. 

Therefore, in order to pick the proper solution of the loop equations one need 
to understand the missing information in the loop equations. When dealing with 
convergent matrix models, the answer is quite simple: the missing information is about 
the choice of the contour of integration of the model. Indeed, we have assumed here 
that the eigenvalues of the matrices were real, but one could take a different one- 
dimensional contour and get the same loop equations as the ones we have derived here 
(as soon as there is no hard edges, i.e. no boundary terms in integration by parts). 
Nevertheless, different integration paths would lead to different correlation functions 
and thus to different solutions of the loop equations. In the case of formal matrix 
models, the path of integration and the associated convergence is pointless and some 
other considerations have to fix the choice of the solution of the loop equations (in the 
hermitian case it is the notion of "filling fractions"). In this article we will restrict 
ourselves (specifically at equation (16. 5p ) to a specific kind of solutions of the loop 
equations because it corresponds to the less technical possible case. 

Let us just mention that, in the hermitian case, the homology space of possible inte- 
gration paths, has exactly the same dimension as the number of unknown coefficients 
of the polynomial PQ°\x,y), and basically, every choice of Pq(x, y) is acceptable, 
and corresponds to a given integration path for the eigenvalues. However, finding the 
integration path from the knowledge of P^ (x, y) is not so easy in general. 

In the non-hermitian case /3 ^ 1, it is not even clear what the homology space of 
possible integration paths is, and a good understanding of the missing information, 
its connections to the initial model and of the generalization of the notion of "filling 
fractions" for our case is still missing and is postponed for a future article. 

6 The Bethe ansatz 
6.1 WKB approximation 

Here, we work at large /3 and fixed N (working at large (3 or small g s when N is fixed 
is completely equivalent since they are related by g s = ttts), so that our integral: 




(6.1) 
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where X = diag(a;i, . . . , xn), Y = diag(?/i, . . . , y^), with the action 



A[X,Y] = I(TrV 1 (X)+Tr^ 2 (F))-2^1nA(X)-2^1nA(F) 
-lxmg 2 s ki{Ip{X,Y)) + 0(g s ) 



(6.2) 



can be computed by standard saddle point approximation, and in particular, we find 
that 



lim Wi (x) 

(3— >oo , N fixed 



N 1 



1=1 



(6.3) 



i.e. it is a rational fraction with N poles which are the saddle points for X, i.e. the 
points such that 



_d_ 

dxi 



A p=oo [X,Y] = 



d_ 

dyi 



A fi=oo [X,Y]=0. 



(6.4) 



Vi=Vi 



Therefore, Wy{x) = % must be a rational function with N poles, i.e. ip(x) is a 
polynomial of degree N. This is the case we will develop below and that corresponds 
to the generalization of the solution presented in [16] in the case of a Schrodinger 
equation. 

However, we would like to stress again that this is not the only possible regime. 
Indeed, we could also try to study directly the "topological regime" where iV — > oo, 



N 



1 - 



remains finite. This regime is 



T fixed, and /? — > in such a way that h 

more interesting for applications to string theory and AGT conjecture 
In this case, ip(x) would not be a polynomial, instead, the coefficients of E(x,y), i.e. 
the choice of W^\x) would be dictated by asymptotic behaviors in accordance with 
the regime studied. This is more challenging, we have done it for the 1-matrix model 
case [T8l[T9j . and we plan to do it for multi-matrix models in the near future. 

Therefore in the rest of the article, we restrict ourselves to the "Dirac" 
case (the terminology "Dirac" is used here to remind that the density of 
eigenvalues is a discrete measure, i.e. a sum of ^-functions, i.e. a "Dirac 
comb"), where N is fixed and — > oo. In particular it correspond to assume 
that ip(x) is a polynomial (whose degree is the number N of eigenvalues). 



6.2 Introduction of the ansatz in the equations 

So, let us look for a solution where ?p(x) is polynomial in x. Then it is clear from 
the recursion relation 1 15. 10p that all the ^(a;)'s are also polynomials. We will use the 
following notation: 



(6.5) 
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So that: 



Y(x) = V[(x) 



T" 1 



N ^ x- Si 

i=l 



(6.6) 



Here we have simply assumed that ip(x) is a monic (remember that ip(x) is only deter- 
mined up to a multiplicative constant) polynomial of degree 7V and labeled (sj)j = i v .. ) jv 
its complex zeros. Note that this ansatz is very restrictive, since usually a linear ODE of 
the kind ( I5.15P does not admit polynomial solution if the coefficients E^(x) are generic. 
Then, from the definition of Uq(x, y), it is also trivial that Uq°1 (x) can only have simple 
poles at the Sj's. Taking into account the term in V^iy) — x in the definition we have: 



A' 



V0<A;<d 2 : <2(*) = E 



8=1 



U k ,i 
X 5 7 ; 



tk + Sk,ox 



(6.7) 



where the Uk,i are at the moment unknown coefficients. 

6.3 Computation of the s^'s and of the u^/s 

Putting this ansatz back into the recursion relations (I5.10p for the U^l(x) and identi- 
fying the coefficient in —^-7, we get for all k and i: (Note that the double poles cancel) 



ufc-i,* - {v[{ Si ) -^E -4— )"w + S ( E " ffe + ) = ( 6 - 8 ) 

In order to solve this recursion we introduce the following matrices of size N x N: 









Bi4 




S = diag(si, . 


■,s N ) 


H 










T 



(6.9) 



and define the vectors : 





/A 




( Mfe.l^ 






1 




Wfc,2 




e = 




, u k = 












\Uk,N/ 





(6.10) 



Then the previous set of equations turns into: 



T ~ 

Uk-x - Bu k = —(t k - 5 k flS)e 



(6.11) 
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Remember again that we have some extra knowledge for k = d 2 and k = d 2 — 1 (Cf. 
(ESD and ([EZD ): 

T ~ 

^ 2 = , u d2 -i = j-t d2 e. (6.12) 

Starting from Ud 2 -i — jftd 2 e and using the recursion relation, we have: 

T ~ T 
u d2 -2 = Bu d2 -\ + —t d2 -ie = — (t d2 Be + t d2 -ie) 

T ~ T 

u d2 - 3 = Bu d2 -2 + —t d2 - 2 e= — (t d2 B 2 e + t d2 -iBe + t d2 - 2 e) 



down to 



u = Bu x + ^rfie = ^ (td 2 B d2 ~ 1 e + i d2 -iB d *- 2 e + ■■■ + he) 

T ~ 

-Bu = —(t - S)e 
In fact all the previous relations can be summarized into the matrix form: 

rp di — k — \ 



or else in components: 



= n tk+ P +iB p e (6.13) 

u ^=nJ2 E Wi(£ P )m (6-14) 

j=l p=0 



Eventually using the fact that —Bu = jf(t Q — S)ewe eventually get: 

M(B)-S)g=0 (6.15) 
On the other hand, from the definition of B we have: 

= e t (Vl(S)-B) , [ B ,S] = ^(ee t -ld) (6.16) 
Therefore we obtain a system of equations determining the roots (sj)i< 



i<N- 



(V 2 \B)-S)e = , e t (Vl(S)-B) = , [S, B] = ^(Id - e"e*; 



(6.17) 

The last system of equations (16.171) gives us a complete set of equations to compute 
all the Sj's. In fact, this corresponds to a Bethe ansatz [5] found in [16] for the arbitrary- 
(3 one-matrix model (where the authors got in the same context: V'(sj) = ; ). 

Remark: (16. 17ft is a set of algebraic equations which determine the Sj's, and the 
solution is in general not unique. The choice of a solution, is also a choice of a saddle 
point for computing our eigenvalue integral at large (3, and thus it is related to a choice 
of integration path for the eigenvalue integral. The number of solutions of (I6.17p . 
coincides with the dimension of the homology space of all possible integration paths. 
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6.4 Rewriting UQ°\x,y) in terms of S and B 

With the previous results, we have that: 

d 2 N 



x 



k=o i=i x Sl 

d 2 AT (fc-fc-l 



^) + § E E E '^^(n^ 



(6.18) 



fc=0 i,j=l p=0 



Moreover we have: 



d 2 d 2 —k—l d 2 k— 1 

fc=0 p=0 fc=0 p=0 



jg(y) - Vj{B) 
y-B 



(6.19) 



so that: 



Uo°\x,y) = x 



Hence in the context of the Bethe ansatz, we can express easily the function Uo(x,y) 
by: 



Ul?\x,y) 



T 1 Vj{y)-Vi{B) 
N x-S y-B 



e + x-V^(y) 



(6.21) 



6.5 Rewriting Pq \x j ij) in terms of S and B 

We start with equation ( 15 .21) giving that: 

y _ V l( x ) + W(yX ) + L d \ u$>\ Xi y) = (y((x) - y){Vi{y) - x) - P (0) (x, y) + T + ^ 



(6.22) 



We can compute the l.h.s. with the help of the definition (16. 7ft : 

y _ V ; {x) + W(x) + ^d x ) Ui°\x,y) = (V{(x) - y){V{{y) -*) + ^ 

d 2 N d 2 N 

Hv - vm) E E - ^ E E 



rp d 2 N 

+ /vEE 



fc=0 



' x — Si TV ^-^ ^-f (x — Sj) 2 

i=l * fc=0 i=l v ' 



Uk,: 



f^tt (x-8i)(x- 8j y 



T N 1 

t/ + =(-^))E- 



rp d 2 N d 2 N T// / x 

WW - rtWM - x) + £ + E E ^ - E E -^y* 



k=0 t=l 
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fc=0 i=l 



iY 



iV ^ a; — Sj 



(6.23) 



Note that we have used the fact that u^ 2t i = 0. With the help of the recursion relation 
for the Uk/s ( 16. 8ft we find: 



2/ _y i '( x) + I y(x) + ^^)f/ ( S 0) (x 



T 

N 



k=0 i=l 
rp d 2 N 



X Si 



fc=0 j^i 



(X - Sj)(x - Sj) 



2/ 



Observe now that we have the identity: 

N 

Uk,i + u k j 



V0<k<d 2 : - V- 



a* Sj) (Si 5j 



TV 

£ 



(a; - Si) (a; - Sj) 



Moreover we have also (using (16.191) ): 

k=0 i=l * 

so that we are left with: 



T V{(x)-V{(S)V>(y)-V>(B) 
N x-S y-B 




(6.24) 



(6.25) 



(6.26) 



(6.27) 



Eventually, we can summarize the results in the following way. With ( 16. 17ft we can 
compute all the s^'s and then using the relations (I6.13P we can get every u^a (0 < k < d 2 
and 1 < % < N). Hence we have found an algorithm to compute explicitly all [7^ (x)'s 
or equivalently every ^(a^'s (and ip(x)), that is to say that we have the solution of 
the loop equation at the dominant order in the studied limit. In particular, assuming 
that ip(x) is polynomial in x is sufficient to determine it completely which indicates 
that the structure of the loop equations are very rigid and hides integrable structure. 



6.6 Yang- Yang Variational approach of the Bethe ansatz 

The Bethe ansatz equation (16.17P can also be obtained from a variational approach 
with a Yang- Yang action. Consider the following functional : 

^S(S,S,A,u) = tvVt(S)+ txV 2 (S)- tr (S AS A" 1 ) -?-\n(A(S)) 
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~ In(A(5)) + I In det(A) - ^ (Ae - e) (6.28) 

where 5" and S* are diagonal matrices of size N x N, A is a, N x N invertible matrix, u 
is an iV- dimensional vector, and remember that e* = (1, . . . , 1). We now look for the 
extremal values of S. First taking the derivative relatively to u gives: 

dS 

— = 0^Ae = e^A- 1 e = e (6.29) 
ou 

Performing the derivative relatively to S = diag(si, . . . , sn) yields: 

|^ = o v;(s t ) - (asa- 1 ) -^J2^— = ° ( 6 - 3 °) 

OSi V / i,i N ' Sj — So 



Performing the derivative relatively to S = diag(si, . . . , s n ) is the same: 

(^4e^=° ( 6 - 31 ) 



8S 

— = & K^) 

dsi N ^ Si - Si 



Performing the derivative relatively to A = (Ajj)^ is more difficult, consider an in- 
finitesimal variation A — > A + 5A, that yields: 

5S - tr (S5ASA- 1 ) + tr (SASA^SAA- 1 ) - tr (AM) + ^it*cL4e 



- tr ((c^A" 1 (jASA" 1 , 5] + ^(Jrf - Aew 4 )^ 
(6.32) 



thus since A is invertible and Ae = e: 

(95 ~ T 

_ = 0^[A5A- 1 , 1 S] + — (Jd-eu t ) = (6.33) 

But since 5" and S are diagonal we have [ASA -1 , S]ij = 0, and the diagonal term of 
(I6.33P gives Vz : 1 = (eu')j.j i.e. we have Vz : Ui = 1 that is to say 

u — e. 

Then a non diagonal element of the relation ( 16. 33ft gives: 

Vz ^ j : (AS A- 1 ) = T (6.34) 

Combining ( 16 .30 j) and f!6.34|) gives that at the extremum we must have ASA -1 = 
B. Then, we can rewrite (I6.32H as: 

5S = - tr (S5ASA~ r ) + tr (SASA^dAA- 1 ) - tr (A5A) + ^u f 5Ae 
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-tr (A~ 1 (6A) ([S^^SA] + L(id- eu'A^J) 

(6.35) 

Therefore the extremum equation |^ — gives also that 

[S,A~ 1 SA} + ^-(Id-eu t A) = (6.36) 

Again, since S and S are diagonal, we have [S, A~ l SA]i ) i = 0, which implies that 
V« : {eu t A) iyi = 1. Since at the extremum u = e, the last equation gives us that at the 
extremum: 

N 

J^A feij = 1 e l A = e* (6.37) 
fc=i 

Let's now introduce the matrix B = A~ X SA. Then taking a non-diagonal element of 
(JE36J, we find: 

VjV i: hi = (A- l SA) itj = N{ ^_ (6 .3 8) 
We can now rewrite (16.301) as: 

Be = A~ 1 SAe = K>'(S)e 
ABe = AV^e 

Se = y 2 '(5)e 

(6.39) 

Therefore if we combine all previous results we find that an extremum 
of S must satisfies the following identities: 

• u = e 

• Ae = e and e t A = e* 

• [ASA' 1 , S] = %(-Id + ee *) and [S, A^SA] + ^{Id - ee*) = 

• ASA" 1 = 5 and Se = V{{B)e 

We recover all the equations of (16.171) giving us a variational way of deriving all quan- 
tities. Note also that by using (I6.30P and (I6.34[) we have: B = A~ X SA leading to 
Se = V((B)e where remember that B = A~ l SA. The last relations show that the 
system is completely symmetric in tilde quantities and non-tilde quantities. Eventu- 
ally the relation Se = V{(B)e allows to determine directly the Sj similarly to the fact 
that Se = V£ (B) e determines the Sj at the extremum. The determination of A at the 
extremum is unfortunately more complicated. Indeed, Se — V{{B)e and Se = V{(B)e 
determines the Sj and the Sj at the extremum and then the matrices B and B at the 
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extremum. The matrix A can then be computed by the relations B = A~ 1 SA and 
Ae = e or equivalently by B = A~ l SA and AE = e. This leads to the system of 
equations: 

vij : o = Y, A,k Z A " + a *av;(s j )-§a 

N 

1 = £4,k (6-40) 

fe=i 



Remark 6.1 It looks like the Yang- Yang action should be the (3 — > oo limit of the actual 
action A[X, Y] defined in (|6.2p . Unfortunately, since the angular integral Ip(X,Y) is not 
very well known, this cannot be proved at the moment. But vice versa, that would give us 
some knowledge of the angular integral I p. 

6.7 Computation of W^fa, x') and W^°\x, x', x") via the varia- 
tional approach 

Since W^-y = , and since we know W±°\x) = "where the Sj's are 

obtained by the Bethe ansatz i.e. from the extremization of the Yang- Yang action S, 
we can easily find wff^ and W^°\ We will see later [7]another way of getting formulas for 
these functions which can be extended by recursion but having the main disadvantage 
that the symmetry of the functions is not obvious at all. 

6.7.1 Computation of W^\x,x') via the variational approach 

Starting from 



= 1 



we get: 



It remains to compute the quantities g^hn, ■ For that we use that the s^'s are the 
extrema of the functional S(s, s, A, u). 

We introduce the variable R = (s,s,A,u) as a "global" variable to avoid detailing 
each of the cases. It is a vector of dimension N + N + iV 2 + N = 3N + iV 2 . Therefore, 
the functional S(R) is a functional of R whose extremum -R ex tr gives us the solution of 
our model. Thus we have for every variation S and every component v. 



= 5 



dS 

dRi 



\R=R 



cxtr . 
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dS\ 3 ^V d 2 S 



s m).. . + E 5S557 .. _ f «« 



oxt r 



(6.41) 

The first equality comes from the fact that -§§:{R = R cx tr) = since we place ourselves 
at the extremum. Let's introduce the Hessian of the functional at the extremum point: 



It is a 3N + N 2 x 3N + N 2 matrix whose shape is given in appendix |X] The last identity 
gives us that: 



d 2 s 



.dS 

I 

dRiJ \ R 



\R—R. C 



d 2 s 



S((R 



cxtr ) i j 



(6.43) 



3N+N 2 

j = l v 1 J ' \R=Rextr 



3N+N 2 



i=i v J 7 



|-R=i?c 



Let's now specialize the last formula for the case when the variation S is 



a 



avi(x)- 



Since the potential V\ only appears in the term Tr V\ (S) in the functional S it gives: 



dS 



dVAx) 



Tr 



x — S 



(6.44) 



Since it only depends on S (and not S, T or «), we see that in the last sum of ( 16. 43ft . j 
only varies from 1 to N, corresponding to the variables s/s (i.e. dVi(x) ^ = ~ (x-s) 2 
and Vj > N : dRj d y i[x) S = 0). Therefore we find (taking 1 < i < N): 




(6.45) 

Note that this last result is only a small fraction of the results contained in the formula: 

3N+N 2 

S((R ext r)i) = - V (H-\J6^) (6.46) 



e («-% OS 

.7 = 1 V J 7 



\R=Ro 



depending on the choice of variation 5 and the component % you take. With the help 
of formula ( 16.451) . it is straightforward to compute W2°\x,x'): 




(6.47) 
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In particular, in this formalism, it is clear that {x,x') ^ W^> (x, x') is a symmetric 
function since "H and its inverse are symmetric matrices. 

6.7.2 Computation of (x±, X2, £3) via the variational approach 

We use the same method: 

Wt\x 1 y 1 z) = ^- ) wi°\x 1 y) 
For any variation 5 and every component k, we have: 



/ \ 3N+N- o ly-i) 



fc=i 

3AT+AT 2 



k=l k 



(6.48) 



We choose 5 = dVl ^ x ,^ and thus 



i.e. 



dRi dRj dV^x") dRi dRj x" - S 



iii<N,j<N 



5U lJ = { (*"-^) 3 — ( 649 ) 
otherwise, i.e. i > N or j > N. 



Moreover, we have 



and 



dm- 1 ) 3N+N2 r) 3 1 

dR k ^ y h >p dR p dR„dR k K 



dR k d*s 



E/ )k,l f)\A (V' _ 



dV^x") v dRidVx{x") Jk ^{x"-si) 2 

Then we easily find: 



KK 3 j — ^ ^ ( x _ s .)3( x /_ Sj )2( x //_ Sfc )2 -T jy ( x _ Si )2( x /_ Si )3( x //_ Sfc )2 

i,j,k=l i,j,k=l 

N (u-A (n-A N 

I 2T y-> v A:.fcV /fc.j _T Cj.j.fc 

JV ^ (i-s i ) 2 (i;'-s ; ,-) 2 (a;"-Sfe) 3 JV ^ (x-s 4 ) 2 <>'-Sj) 2 (z''-Sfc) 2 

i,j,k=l i,j,k=l 

(6.50) 
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where the C^-^'s are linked with the third derivative of the functional S: 

37V+iV 2 „ 3 _ 

c.,, d = r £ MeM^ifflHr (6 - 51) 

It is again clear that (x,x',x") is a symmetric function of its 

variables. 



7 The topological recursion 

In the last sections we have seen how to compute Uq (x, y), PQ°\x,y) and W\ (x) in 
the context of the Bethe ansatz. Then by a recursive application of the loop insertion 
operators d y x ^ we can compute easily all the Wn s Un's and Pn s. Thus, for now 
we have all desired quantities at the dominant order (i.e. g s = 0) but we still miss the 
subleading corrections g > 1. It is the purpose of this section to present a recursive 
algorithm to get all the subleading corrections. 

The algorithm is very similar to the "topological recursion" considered in [4TRI45P48] , 
and is in fact more comparable to another version of the topological recursion for the 
hermitian 2-matrix model, described in 



7.1 Higher loop equations 

Evaluating the coefficients of g 2 s 9 in (I3.28P gives the higher loop equations for g > 1: 



(y - Y(x) + ^d x ) U^(x,y) + w[ 9 \x)uf\x, y) = -£w^\x)U< h \x,v) 
-Ut^ix, y; x) + pxUt^ix, y) - P (s) (x, y) - f 5 g ^ 



(7.1) 

We can also project the higher order loop equations (13.391) onto g 2 s 9 (with n > 1) 
to find (we note £ = {£i, . . . , £ n } for short): 



y - Y(x) + ^d x \ U^(x,y;i) + U<£\x,y-§w[ 9 \x) 
-■-J2ui h \x,y;£)Wt h \x) + ^d x U^(x,y;i) - (x,y;x,g) 

h=l 

9 nr. 

E E Wil- I ]\xj)ui% l (x,y;t\l)-Pi 9 \x,y;0-J2- 
h =°ic£,i^9 k= 

U^_l(x, y, £1, ■ ■ ■ , £fc-i, £fc+l; • • • , £n) — Un-l(£k, V\ &b • • • ; £fc-l; Cfc+b ■■■i^ 

x-ik 



t 9 ^ 



(7.2) 
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We shall also write those equations for each power of y , using 

d 2 -l 



k=0 
d 2 -l 



(7.3) 



(7.4) 



fc=0 



7.2 Matrix form of the loop equations for g > 1 

The goal of this subsection is to rewrite the loop equations (17. ip (0 < k < d 2 + 1) 



T 



-Y{x) + jjd m \U$(x) 
U[% 1] {x- x) + ^d x U^ l \x) - ^6 g=1 6 k -_ 



5-1 



(7.5) 



-J2 W i 9 - h \x)U$(x)-P$(x) 



h=0 



into a matrix form suitable for our algorithm. 

First by taking the case k = d 2 + 1 and k = d 2 we obtain that: 

\fg>0: U$ 2 (x) = , U^ix) = t d2 Wi 9 \x) 

Now let's introduce the following d 2 x d 2 matrix: 



/ 



V(x) 



-Y(x) 
1 



... 
-Y(x) ... 



\ 



V 



1 -Y[x) 



U {0) (x) 

*d 2 
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Then we can rewrite the loop equations (17. 5p into the following matrix form: 



d x + V{x) 



f- 1 rr(Q ( t ) 



1=1 



t 



\U^(x)J 

(u!> 9 l) (x)\ 
U^ l \x) 



( P${*) \ 



( U[ 9 o 1] (x;x)\ 
u{ 9 i 1] {x;x) 



\u[ 9 d2 1 } 1 (x;x)J 

fl\ 



w 



(7.6) 



(7.7) 



(7- 
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In the same way, we can also rewrite the projection of the higher loop equations 
17.21 onto y k as: 



vSH-ite 3 + (^9 X - Y(x)) U®(x; |) = - £ <fe O^i 3 "^' 

^ ' ft=0 



j=l M 

^n-l.fcO^ Cl? • • • ) • • • ) Cn) — ^n-l,fc(£j> £1, . . . , £j+l> • • • ) 



>(<?), 



(7.9) 



Evaluating the coefficients in k = + 1 and k = c? 2 leads to: 



v« > i , > o : u { n 9 Ux,£) = o , ulfX^o = uMU^O 



(7.10) 
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and thus leads to the linear system: 



^d x + V{x) 



^n+l,l( X ! X i 



£ 



h=o /uJ r =£/(/,/o?H(0,fl)u(So)} 



( ^i\ [0 ( a ,gl,...^-l^ + l,...,g„)-^ 9 _ ) 1[0 (g,-,6,---.€j-l.ej + l,---,gn) \ 
^i\a(^.--^-l^ + l.--,gn)-t/^ ) 1 , 1 «^6v--^-l,gj + l,---,gn) 



(g) 



-1,^2- 



1 (x,6,-,g J --l,g J - + l,-,gn)-^ 9 _ ) 1[d2 _ 1 fe,6,-,^-l,gj + l,-,^) 



(7-11] 



Observe now that the r.h.s. of the last two systems ( 17. 81) ( 17. lip is composed of 
lower orders in g and an unknown polynomial vector Pjfl- The l.h.s. is always the 
same at every order in g and for every n. Therefore, if we can find a way to get 
rid of the unknown polynomials P^l's in the inversion of the system, we see that we 
can perform a recursion and get all the corrections recursively. We will do so in the 
following sections by introducing a suitable kernel K(xq,x) and using some residue 
methods (which will automatically get rid of the polynomials). 

7.3 Inversion of the linear system: kernels K(xq, x) and G(xq, x) 
Suppose that we can find matrices K(xq,x) and G(x ,x) of size d 2 x d 2 such that: 
• They satisfy the relation: 



(v\x) - jjd^J K(x ,x) = G(x ,x) 



(7.12) 
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• The function G(xq, x) has the form: 

G(*o,*) = ^ + y^) (7-13) 

Xq — X X — Sj 

j=l J 

where the A^a^'s are d<i x <i 2 matrices. 

• The function x i— > K(xq,x) is analytic at every x = Sj 

If we can find such matrices, then we claim that we can invert any system of the 
form (I7.8p : (jfd x + T>(x)) u(x) = v(x), where v(x) is assumed to be known up to a 
polynomial and u(x) is assumed to have only poles at the Sj's and behaves like O (^) 
when x — > oo. Indeed under these assumptions we have: 



— Res u t (x)G(xo,x 

X—>X() 

N 

ERes u t (x)G(x , 
X — >Si 

i=l 

N , 

V Res u\x) V t 

i=i v 

jr Res u\x) (v\x) + ^-t) K( 

i=i x s% V 

TV 

ERes w*(a;)-ft'(a;o, x) 
X^Si 



T 

x) - —d x ) K(x , x) 



Xq,X) 



(7.14) 



Note that in the last equality, the polynomial part of v(x) does not contribute only 
the terms in v(x) which have poles at the Sj's contribute. Hence we have successfully 
inverted our system: 



N 



U 



x ) = y~] Res if'Oo, 

i=l 



x)v(x) 



(7.15) 
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7.4 The topological recursion 

Therefore with the help of the results in the last section we can compute all correlation 
functions by the recursion: 



( U[% l \x-x)\ 
u{ 9 ^\x-x) 



9-1 rr(0 



N g-l 

i=i 1 i=i 



KX-iW 



+ 



( U& l \x)\ 

[y^x)/ 



(9-1) 
l,d 2 -l^' 



x: x 



N 

Y 



Or 



(U^ 1] {x)\ 
U^ l \x) 



V 



(7.16) 



or more generally for higher order correlation functions: 



U^(x ,0 



V Res KHx q , 



( U^{x-x,Q \ 

^n+l,l( X ' X i 



i=l 



(U^ix;^ 



V^" 2 -i(^oy 



E E 

^=o juj=£/ (J,fc)7q(0, fl )u(£o)} 



(ff-fc). 



x, J) 



^ U^lifiix, £l> • • • >0'-l'0+l> • • • >£n) ^ 
^n-l,l( X ' 0> • • • ) 0-1' 0+1' •••)&*) 



V^n-M 2 -l( X ' 0' • • • ' 0-1' 0+1' ■■■■>£n)J 
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7.5 Determination of the kernel K(xo, x) 

The determination of the kernels K(xq,x) and G(xo,x) is presented in appendix [Bj 
We present in details some recursive formulas to obtain the various K(xq, Sj) and the 
derivatives K'(xo, s»), . . . , K^ n '(xo, si) in the appendix. The results are the following. 
If we define the following vectors (of size e^): 



o 

V 1 / 



/ 



W; 



U , 



(7.18) 



\ u d 2 -l,i / 

Then the matrices A^xq) defining G(x , x) are connected to the K(x , Sj) by: 




The matrices K(xq, Sj) are determined as a solution of the following system: 



MK q = f q 



where 



/ \ 



K„ 



( K(x ,si) , q \ 
K(x , s N ) 0tq 
if (z , si) I,? 
K{x , s N ) l q 



K{X Q , Si)d2-l,q 
\K(x ,S N )da-l,qJ 



) Jq 



XQ—Sl 



xq—sm 







(x'u-si) 2 

I u q,N . 

\ (x -s N ) 2 ' 



( o \ 



^2 



1 



\ (xo-s N ) 2 ) 



(7.19) 



(7.20) 



(7.21) 
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where the Xf> 1 _ s in f q are in position Nq + i. (Note that fd 2 -i does not have these 
terms). 



M 



( —B* Id 
-B l Id 





\ 





(7.22) 



-B* Id 
\ B Bi ... B d2 _ 2 Bd 2 -iJ 

where all the matrices Bk, Id and are symmetric matrices of size N x N given by: 



Bo= N 



( (£0)1,1 

"0,l+ttQ,2 
(si-s 2 ) 2 



"0,iV+"0,l 
\ (SN-Sl) 2 



"0,1+MQ,2 

(si-s 2 ) 2 



U 0,N+ U 0,N-1 
(S]v-sjv-l) 2 



UQ,1+UQ,N \ 

(si-sjv) 2 



U 0,N-1+ U 0,N 

(s]v_i-S]v) 2 
(Bq)n,N J 



(7.23) 



with (-Bq)m = - 
And VA; > 0: 



B, 



T K 1 l S iJ M 0,i - 2^ (s,-s,) 2 



1. 



I 1-DfcJl.l ]V ( Sl -s 2 ) 2 



JV (si-s 2 ) 2 



T "fc,jv+"fc,i 

\jV (SAT-Si) 2 



with = -W(aO«fc,i - JjE^f- 

37* 



T "fc.iV+Mfc.JV-l 

iV (sjv-sjv-i) 2 



r ^fc,i+"fc,jv 

N (si-sjv) 2 



T "fc,JV-i+Mfc,jy 
TV (s]v_i-sjv) 2 

(Bk)N,N J 



(7.24) 



If we define 



d 2 -l 

j=k 



we have 



/ o 
Id 

B f 



0\ 



B 



Id '•• 

B l Id 



\5* d2 " 2 ... 5 f2 B l Id 0/ 



/ Id \ 

B t2 



\B td2 ~ l ) 



C^id C 2 ... CVi W) 



(7.25) 



That gives if q < d% — 1: 

JW*o,*) = E 0(p > q) 



X - Sj 
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J2(B"C^C, kj -L- - 2(B»cr') u7 -^j (7.26) 



Xo - Sj J (x - Sj)' 

j 3 

where 9{p > q) is the heavyside function, equal to 1 if p > q and otherwise. And 
when q = d 2 — 1: 

JiWi(so, s «) = ** E^o" 1 )^ , 1 (7.27) 

iV (^Xo - Sj) 

Eventually the higher derivatives in x of K(xq,x) evaluated at x = s, are obtained 
recursively with formulas flfOOj) . flR32j) . flR33j) . flR39j) and f lR40j) derived in details 
in appendix. Note also that to complete these formulas, we need to prove that the 
system of equations is consistent. This involves a special identity which is proved 
in appendix O As a conclusion, we have obtained here a recursion to com- 
pute all correlations functions Un\x,y; x). This corresponds to a generalization 
of the topological recursion jlQllIE] for the arbitrary-/? two-matrix model studied here. 
Moreover, when d 2 = 2, one can recover the case of the arbitrary-/? one-matrix model 
developped in [46] . 



7.6 Leading free energy 

We clearly have from the saddle-point approximation method, that 



/o = -^5. (7.28) 



where the Yang- Yang action S 



S = trV 1 (S)+ trV 2 (S)- ti (S AS A' 1 ) -^\n(A(S)) 

~ HA(S)) + ^ In det(A) - ^tf (Ae- e) (7.29) 

is evaluated at its extremum. 

The link H4.8[) between the WKB expansion and the standard topological expansion 
gives 

F (^) = fo(N) = -S (7.30) 

7.7 Subleading g = 1 

To the first subleading order, our topological recursion gives 

= V Res K\x ,x) (u[°\x;x) - ^d x U^(x)) (7.31) 

i ' 

and we remind that the d 2 — 1 th component of the vector (xq) is the function 
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In order to perform the computation, we first need to compute u[°\x; x), which is 
also computed by the topological recursion and is worth 

N 



^Res K\x Q ,x)Ul?\x) 



(7.32) 



i.e. 



N d2_1 N 1 

i=0 i=l ^ 1 ^' 



Lk UlA 



(7.33) 



and we need it at x = £ i.e. 



d 2 -l N 



:Ul, 



(7.34) 



2=0 i=l v 1 ' 

Notice that K(x,Si) may have double poles at x — Sj, and thus U{°l(x,x) may have 
poles of order up to 4. This implies that (xq) can be expressed in terms of K(xq, Si), 
K'(xo, Si ), K'"(x , Si ). 

From the WKB approximation, we know that we must have 



/ 1 = -_lndet.A" 



(7.35) 



where A is the action (see (I6.2p ). which is approximated by the Yang- Yang action S 
in the limit f3 —> oo, and Ji is the Hessian of the Yang- Yang action S, so this leads us 
to conjecture that 

(7.36) 



/i = -~lndetH. 



We find that: 



Wi 1 \x) = 1 ^ 

tdo 



B, 



(X 



(7.37) 

where we could explicitly verify the lack of a simple pole term. The coefficients Bj are 
given by: 



i.n 



p=l q,n,m x v ' 

+2 % V- 1 ]* 



(7.38) 



The technical details (and the definitions of quantities involved in the last formula) of 
this computation can be found in appendix [E] with other possible expressions for Bj 
given in flK68|) or flK68|) . 
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7.8 Example 

As an example, we present explicit computations of the method in the case when N = 1 
in appendix (jD]). In particular in this case we are able to determine W^\x) explicitly 
and to compute the corresponding F\. In particular we show in this case that it satisfies 
the equation: 

/ 1 = --lndet7f (7.39) 

7.9 WKB computation 

The WKB approximation, consists in expanding the integrand in the vicinity of a saddle 
point. To leading order, one gets a Gaussian integral, and corrections, are moments 
of the Gaussian integral, and they are computed by Wick's theorem. In principle, one 
can write: 

Z= f dXdYe~^ A[X,Y] (7.40) 



and expand: 

X = X + g s X , Y = Y + g s y (7.41) 

and expand 

A[X, Y] = A[X, Y] + l -gl A"[X, y] + gl 6A[X, y] (7.42) 

where the Hessian A" is the quadratic form of second derivatives of lim A, and 5 A can 

be computed by its Taylor expansion in g s . 
We thus get 

Z = fiXi^i™ ~ g^i^ (7.43) 

J VtetA 1 $ dxdy 

i.e. we recover 

/o = - lim A[X, Y] (7.44) 

p— >oo 

and 

fi = - lim - det A", (7.45) 

0— >oo 2 

and in principle, a standard WKB approximation using Wick's theorem, should allow 
one to find / 2 , / 3 , . . . explicitly. Then, wi g) can be recovered by applying repeatedly. 
However, this is long and tedious, and difficult to do in a systematic way. One of the 
difficulties, is that one has to compute the large (3 expansion of Ip(X, Y), which is not 
very well known at the present time. 

So, our loop equation method is an alternative to WKB, it can be performed in 
a systematic way, and gives all orders in the g s expansion. Also, notice that Wick's 
theorem gives an expansion of Z in powers of g s , whereas our method gives directly 
the expansion of In Z. 
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8 Conclusion 

We have presented here the loop equations of the two-matrix arbitrary-/? matrix models 
and a recursive way to solve them for a specific polynomial Bethe ansatz. This work is 
a generalization of the arbitrary-/? one matrix model with the same kind of ansatz as 
developed in jl6]. This work is also part of the "quantum algebraic geometry" project 
started for the one-matrix model in [T8jll9p46j . In particular, we have proved that in the 
two-matrix arbitrary-/? matrix model, we find a "quantum curve" (15.151) of arbitrary 
degree. For now we have only solved this model under the assumption of a polynomial 
ansatz for the solution corresponding to a specific limit of the model (N and T fixed 
while g s — >■ (i.e. /? — > oo)). However, we expect the notions developed in [T8|[T9] 
to generalize well for the two- matrix arbitrary-/? matrix model as it has been the case 
here. This future work considering generic solutions of the "quantum curve" (I5.15P is 
under progress at the moment. Then one may wonder if properties of integrability and 
of spectral invariants [HH2D1SS] known for the hermitian case can be generalized in the 
arbitrary-/? case. 
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A Hessian matrix of the Yang— Yang function 

Using the notation of the paragraph regarding the variational approach and using 
(I6.30p . we have: 




(A.l) 




(A.2) 




(A.3) 




(A.4) 
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/ d 2 S 
Vl<i,r,s<N: U^- 



N 



|extr 



$r,i s s(A 1 ) s ,i ~ Aj iP S p (A 1 ) PjT .(y4 1 ) 



p=l 

A? 



p=l 



(A.5) 



Where we have used here: 



(A.6) 



We can do the same kind of computation starting from (16.311) : 

d 2 S 



5.^2 



Vl< % < N : 



dsids t/ |extr 



VI < i,j < N : 



|extr N & ~ ~ S i 

d 2 S 



dkdujs |cxtr 



(A.7) 

(A.8) 
(A.9) 



VI < z, r, s < A r : 



d 2 s 



dsidA r , 



JV 



|extr 



&s,iSr(A )i,r ^ ^ Ap^Sp^A ^ S p(A , 



p=l 



(A.10) 



The computation relatively to u is easy: 



VI < i,r,s < N : 



«9 2 <S 



duidA r 



s / | extr 



VI < i,j < N : 



d 2 S 



(- 



d 2 s 



( d 2 s 



. du i ds j J |cxtr V duidSj J |extr V^a^ J |extr 
Eventually the most technical one is the double derivative relatively to A: 

d 2 S 



(A.ll) 
(A.12) 



V 1 < a, b, r, s < N : 



8A a hdA r 



extr 



Eventually the shape of the matrix looks like: 



V / b,r 

(A.13) 



ft 



/GnxN Xn x n Y NxN 2 Oatxjv\ 

X Gnxn Z NxN 2 Onxn 
Y l Tn 2 xn 2 In 2 xn 

OjVxJV / 



\0jVxJV ®NxN / 



(A.14) 
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Where the matrix G is given by: 

Gii = 



and the matrix G is its dual: 

Gii 



V{{ Si ) 



T 
N 



y— 



T 



N( Si - Sj y 



(A.15) 



T 



Z-^ ( q . _ S ■ 



Z.\2 



(A.16) 



The matrix In 2 xn is is a A 2 x N matrix with N times — ^ in each col- 
umn and zeros elsewhere. If we sort the A r ^s in the following way: 
(Ai t i, . . . , Ai t N, . . . , A 2> i, ■ ■ ■ A 2; n, ■ ■ ■ , An^n) then the matrix looks like: 

(e 0\ 

e ■■■ 



' N 2 xN 



T 

'N 



(A.17) 



e 

\0 ... e/ 

The other matrices X, Y and Z are expressed by the previous computations, but since 
they do not have a compact or interesting form we do not mention them here. 



B Finding kernels K(xq,x) and G(xq,x) 
B.l Decomposition of the matrix V t (x) 

Before presenting the way to find the kernels K(xq,x) and G(xq,x), we need to in- 
troduce some notations. First, since T> l (x) is going to play an important role, we 
decompose it in the following way: 



V\x) = -Y(z)Id + A + ^rU\x) 

td 2 



where we have defined the vectors and matrices of dimension d 2 by: 

/0 1 . 



M 
o 







/ Ufl(x) \ 



u 



A 



\o ... 



i 

0/ 



(B.l) 



(B.2) 
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Note here that all the matrices defined here are already known since they only 
involve leading order quantities determined earlier. In particular we have: 



T N 1 
Y(x) = Vl( X )--J2 — 



N ^— ' x — Si 

i=i 1 



(B-3) 



and: 



N 



i=l 



X Si 



(B.4) 

where the coefficients u^^ are known and satisfy ( I6.13p . In the following equations, we 
will also need to write down a series expansion of U around x = Sj. Therefore, it is 
natural to introduce the following vectors: 



W; 



1 

td 2 



^ UQ,i \ 



Vi 



1 

td 2 



E 



1*1, 



-U 



E U d 2 -l,j _ ± 
Si -S4 L d 2 -1 



(B.5) 



V} 



1 

td 2 



. ,.(*-«j) 2 



"d 2 -i,j 
V ft'.Csi-aj) 2 



Vr 



1 

*<f 2 



/y "o,j \ 



) 



(Si-Sj) 3 



u d 2 -i,j 



(B.6) 



which correspond to the first terms of the expansion: 



t d2 U 



w, 



X S4 



+ V, + VHx - sA + VHx - s t ) 2 + 0(x- Sl 



(B.7) 



Remember that u d2 -\,i = 4?t d2 so that we have the following important relations: 



(B. 




B.2 Solving the system: Finding G(xq, x) and K(xq, x) 

In order to apply our scheme of recursion we do not need to compute explicitly K(xq, x) 
but only to be able to determine K(xo, Si), K'{xq, s») and so on (where a prime means 
a derivative relatively to the second variable: K'(xq,x) = d x K(xo,x)). Indeed, if we 
know this quantities then we can compute all the residues we need in our scheme and 
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invert the differential system. This is reassuring since it was clear that K(xq, x) was not 
uniquely defined since it is a solution of a differential equation. We will first determine 
the kernel G(x,Xq) by computing the matrices A^xq) presented in (17. 131) . Note that 
in order to have simpler expressions, we will label the components of the 
matrices K(xo,x) and G(xq,x) from to g?2 — 1 instead of the standard 1 to 
d 2 . 

The strategy is the following: we will identify the coefficients of order (x — Si) k in 
equation (17.121) . 



Terms in (x — s,0 



-l. 



A straightforward computation gives: 




Multiplying on the left by wl and using (IB. 81) we find that: 

T 

2—w i t K(x , = w i A i (x Q ) 

so that: 




(B.9) 



(B.10) 



(B.ll) 

Therefore, we see that the knowledge of K(xq, s,) is completely equivalent to the knowl- 
edge of Ai(xo). In terms of components we see that only the last line is special: 



Vg < d 2 , Vp < d 2 - 1 : {A i {x )) m = — (K(x , 8i )) pi 



(B.12) 



d 2 -2 



Vq<d 2 : (Ai(x )) d2 _ l q = — (K(x , Si)) d2 _ lq + ^Y1 u ^ ( K ( x °> s i))k, q ( B - 13 ) 

td * fc=0 

We need to go to the next order in order to determine K(xq, Si) or equivalently 
Ai(x ): 

Terms in (x — s,)°: 



Looking at the constant coefficient of the expansion gives: 



V/ ) K(x , s i )+W t K'(x , Si ) = + A 



Xq — Si 



Si Sj 

(B.14) 
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Terms in (x — si) . 



j# {Si Sj) 



2 ^ Vo , s , ) + w ~.^_ 5 _M^L (B , 5) 



We see that the first equation gives can determine the first lines of K(xq, as soon 
as we know the last line of K(xq, Si). In the last equation, we see that we can get rid of 
the second derivatives by multiplying on the left by fw\ and using wjf= S. It gives: 



jj-i iV \ S i S j) J iV 

The last equation is only interesting on the last line, because all other lines trivially 
give zero (because we multiply on the left by r). Now we note the following identity 
for the component (c^ — 1, q)' 

irw\ (-Vl(si) + J2 Nis ^_ s .) + A j ^ + ^K'(x , si) 

\ \ * / / da-1,9 

-^i(si) + E V E T^ K '{ x v, si) k>q + V -^^'(xo, si) h>q 

N[8i- 8j )J k=Q t d2 k=i t d2 

T d2 ~ l ( u k \ 
+77r~ Yl — — — *k + «i4,o K'(x , sij^ q 

= ( U k-hi ~ K( S i) U k,i + T? E Uk ' 1 _ Uk ' 3 + A/^ - ** + s ^k,o))K'(x 0l Si) k 

td 2 k=Q ^ j¥H Si Sj 1\ 

= (B.17) 

Where we have used the relation for the u k /s (16.81) . Therefore, we can simplify 
(IB. 161) into a much simpler way: 



\-V?{si) - N{ *_ 8 .y ) s ^ + ^'Kixo, s t 
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(x - Sif 



fw t i A j (x ) 



[Si - Sj, 



(B.18) 



This time, we see that all the derivatives K'(xq, Sj) have vanished. Observe that we 
can replace A,(x ) by some K(x , Sj)'s using (IB.9j) . 



fwj -Vf( Si ) - £ 



T \ T 

] if (x , s{) + — r^ u if (x , 5 4 



N{si - Sj ) 



T sr^ + rw t j )K(x , Sj) 



N 



(x - s t ) 2 N 



E 



(B.19) 



Taking the component (d 2 — 1, 9) of the last equation we have: 



fc=0 



It 



TV 

d 2 -l 



3& 



\S{ Sjj 



q.i 



(X - S, 



_T >p (uk,i + u kd )K(x , Sj) k: g 

\2 AT 2-^i 2-^i 



,Si - Sj) 



(B.20) 



k=0 j=£i 

Remember that from ( IB.14I) we have some recursion for the other lines. Vp < d,2 — 1 



-V[{si 



5. 



p.q 



Xq — Sj 



(B.21) 



Therefore, we can group (IB.2Q[) and (IB.2ip into a matrix form. Let's introduce the 
following vectors: 

/ \ 



K„ 



( i^(x ,Sl) ,g \ 

K(x , S N ) 0: q 

K(X , Sl)l ( g 

K(x , s N ) hq 

K(X , Sl)d 2 _l,g 

\K(X , S N ) d2 -l,q/ 



XQ~S N 



(xo-si) 2 



V (x -s N ) 2 
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/ \ 



) fd 2 - 



T ~ 
N 



td 2 




(B.22) 



(x -S N ) 



_ J 



where the are in position Nq + i. (Note that fd 2 -i does not have these terms). 



(-B 1 Id 
-£* Id 



\ 





... -B l Id 

\B Bi ... B d2 _ 2 B d2 _iJ 



(B.23) 



where all the matrices B^, Id and are symmetric matrices of size N x N given by: 

/ ( D \ «0,l+ u 0,2 



#0 



r 

iV 



(si-s 2 ) 2 



"o,i+"o,jv \ 
(si-s N ) 2 



"0,l+"0,2 
(si-S2) 2 



"0, AT +"0,1 



U0,N-1+UQ,N 



U 0,N+ U 0,N-1 
(SJV- 



(B.24) 



Tp- (Bo)n,n J 



with (JB| 



And Vfc > 0: 



T "fc,l+"fc,2 

N (si-S2) 2 



T "fc,JV+"fc,i 
\JV ( SJV - S1 ) 2 



T "fc,l+"fc,2 

JV (si-s 2 ) 2 



T "fc,iv+"fc,jv-i 
N (sjv-sjv-i) 2 



T "fc,i+"fc,JV \ 
N ( Sl - Sjv ) 2 



T "fc,iv-i+"fc,iv 
TV (sjv-i-Sjv) 2 



k)N,N 



(B.25) 



with CBt 



V\ \S % )U k ,t Jjl^ ( Sl - Sj ) 2 



Then, we can rewrite the equations (IB.2Q[) and (IB.2ip into the nice matrix form: 



V0 < q < d 2 - 1 : MK q = f q 

(B.26) 

Hence, by a simple matrix inversion, we can find all the K(xq, si) and then using 
(IB. 9 1) find the kernel G(xq, x). Note in particular that the last column of K(xq, si) and 
hence G(xq,x) only have terms in r Xo _ 8 .\z ■ In particular, if we define the "Bergmann 
kernel" B(xo,x) = —\d x G{xQ,x), we observe that it is a symmetric function. Eventu- 
ally, note that in the case d 2 = 2 (i.e. the one-matrix model) we recover the results of 
the one-matrix model as described in [46] since we have a scalar problem corresponding 
to the component (d 2 — 1> d% — 1) of our problem. 



B.3 Computation/ choice of K'(xq, Sj) and K"(xo : s$) 

Now that the i^(xo,Sj)'s and the kernel G(a?o,x) have been determined, the next step 
is to find the higher derivatives K'(xq, si), K"(x , si), . . . in order to be able to apply 
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the residue formula (17.151) . Unfortunately, since K(xq,x) is defined as as solution of a 
differential equation, it is not determined uniquely. For instance we can write the loop 
equation projected at orders (x — Sj)°, (x — Si) 1 and (x — Sj) 2 . Using notations (IB. 5|) 
we find: 
order 0: 



A - V{( Si ) + L V ^— + rV/ ) K(x , s i )+W t K'(x , s t ) = + T 

N ^ Si - sj / x - Si *rf. Si - Sj 

(B.27) 



order 1: 



order 2: 



-Vfl*) + E s .) + A + ) A"^' s -> 



—if (s , a .) + ™, - — — - - ^ _ — - (B.28) 



+E ^ + * + ^ 

^w-E^^+^'W*) 

3iV 6 (s - Si) 3 ^ (s< - s^) 3 

Remember that if we multiply on the left by fw\ we kill the terms 



>t 



-V{(si) + N ^_ s a + A + rVj J . We clearly see that now that the K(x , Sj)'s are 

known, we have some arbitrariness in solving this set of equations. Moreover, it is not 
obvious that this set of equations is compatible with our determination of K'(xo, Si)'s 
in the last section. Fortunately, we prove in appendix [C] that this is the case. For 
our residue scheme, we can choose whatever kernel if(so,s) we want providing that it 
satisfies the assumption described earlier. Here we will choose it so that: 



Wp < d 2 - 1 , \fq,i : K'(x , Si) p>q = , Wq, i : K"(x , Si)d 2 -i,q = 



(B.30) 
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We claim that this assumption is compatible with our previous set of equations. 
Indeed, since the matrices K'(xq, Sj) have vanishing d 2 — 2 first lines, we can extract 
from (IB.27p . the missing components: 



N 



T 



N{ Si - 8j ) 



K(x 0l Si) d2 - h , 



°q,d2-l 



X - Si 
Sq,d 2 -1 



EE 



Uk,j 



Si Sj 



h + Sk,oSi K(x , s 



i)k,q 



x - Si 
d 2 -l 



k=0 

+ V{(si)K(x , Si) d2 _i i9 



T y^K{x ,Si ) d2 ^ q - K(X , Si),^-!,, 



N 



Si Sj 



+ I X^T^rV ( K ( x Q> Si)k,q ~ K(x , Sj) k>q ) + (-t k + 5 k , Si)K(x , Si) k , q 

(B.31) 



This expression is also equivalent to the following (multiply (1B.27j) on the left by 
fw\ and use ( IB. 171) to get rid of K(xq, Sj)): 



(S) 2 K '( X 0, Si)da-l, 



rf 2 _1 A I \ 



+ EE- 



q t d2 (x -Si) 1 < id 2 (si-Sj) 
k—\i jt=2 

di— 1 

Ug,i _|_ T_ ^2 { U k,i+ U k,j)K{XQ,Sj) ktq 



td 2 { x o-Si) N 



k=0j^i 



t d 2 ( s i~ s j) 



(B.32) 



Then we can compute the first lines of K"(xq, Sj) using ( IB. 281) : Vp < c?2 — 1 : 



2N^ i^Oi s i)p,q 



5 P)d2 - 2 K'(xo,Si) d2 ^ q - V"(Si) + E 



N( Si - Sj y 



K(x ,Si) 



p.q 



(XQ-Si) 2 



I A j (x ) Pt g 

37* 



5 Pyd2 - 2 K ; (x , Si) d2 -. 1>q - V"(si)K(x , Si) p . q - 



T K(x ,Si)p t q-K(xo,Sj) P: q 

N ^ 



(X -Si) 2 



(B.33) 

where the r.h.s. is known. We see also that in (1B.28j) . the last line of K"(xq, Sj) vanishes 
identically so it cannot be determined and we can choose it to be zero to simplify the 
computations. These degree of freedom were also present in the one-matrix model [^6] 
where K"(xo,Si) (which was scalar) remained unknown. Moreover, if we count the 
degrees of freedom, we see that we have exactly d 2 2 unknowns in the determination 
of K(xq,x) which corresponds to a unknown matrix of size d 2 x d 2 as expected for a 
linear differential equation of the form (I7.12p . 
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B.4 Recursion to determine K^ u \xq : s$) 

Now that we have the matrices K(xq, s$), K'(xq, Sj) and K"(xq, Sj), we can perform a 
recursion to get the higher derivatives. To obtain them, one needs to look at the terms 
of the expansion of order {x — Si) n . The main problem here is not conceptual, but 
rather notational. We introduce the following expansions: 



Y(x) 



T 



N(x - Si ) 



p=0 



pi 



N(s t - 



^ ] Y(x-si)p {x ~ s i 

p=-l 

r(0) 



(B.34) 



A 1 



From the definition of Uq fe (x) = ~ tk + 8k,oX we have also: 



X Si 



E 



(-1)% 
( Si - 



E 



+E E 



r(0) 



+ 4,0 (a; - 



P =-i 



We can regroup them into the following vectors: 



(x-Si)P 



( u m \ 

1 U 0,0,(x-Si)P » 
rr(0) 



\ U {0) I 

\ u d2-l,0,(x-Si)P/ 



Note for example that we have in the former notations Wi = U { 
V} = U, 



(B.36) 



(x-Si) 1 y 



We will also need an expansion in (x — Sj) p of the x i— > G(xq, x) function: 

Mx ) id , ^(-i)^^)' 



G(a;o,a;) 



E 



E 



OO 



P =-i 



(B.37) 



Now it is easy to write the equality of the power n > 1 of the differential equation 
{V\x) - jjd x ) K(x Q ,x) = G(x Q ,x): 

Terms in (x — Si) n : 



T K( n +V(x , Si ) ^ KW(x Q , Si ) K^(x , Sl ) 
-n— — ; — — — > J Y {x _ s . )n -i - h A- 



N (ra + 1)! 



1=0 



/! 



??! 
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+m (n+1)! + Vh iX ~ Si) 



l\ 



(B.38) 



We can project this equation on the first — 1 lines. Any term with a r does not 
contribute. We find for p < d 2 — 1- 



n 



iV (n+1)! _ ^(x- Si )"{Xo)p,, 

+ 



(n+1)! ~~ ^^-s i j"V a 'U / /p,g 



1=0 



[B.39) 



Note that the r.h.s is completely known by recursion since it involves derivatives 
only up to order n. Then since we know the first lines, we can compute the missing 
last one by projecting the equation on the last line: 



V 1 l )N (n+1)! ~~ 

KW(x ,s i ) d2 _ 1 . g 

1=0 

di — \n 
'"E~ S Zl^fc,(x-s I )"- i 



-G(x-Si) n [Xo)d2-l,q 



d 2 -2 



+ £ 



',i 2 



k=0 1=0 



k=0 

K {l HxQ,Si)k,q 
l\ 



(n+1)! 



^B.40) 



Again, we note that the r.h.s. is known by recursion since it involves only derivatives 
up to order n. Moreover the coefficient (n — 1)-^ is never zero when n > 1 so that 
we can compute {^o, Si)d 2 ~i,q- Therefore we can determine with both ( IB. 391) 

and ( 1B.40I) all the components of K^ n+1 \xo, si) as soon as we know the lower orders. 
Since we know K(xo,Si), K'(xo,si) and K"(xo,Si) we can compute by recursion 
every derivative K^ n \xQ,si) and thus implement the scheme given by ( 15 . 1 Of) 
to determine all functions U[fl(x) and then by taking the last component of 
this vector (17. 6ft the correlation functions w[ g \x). 



C Compatibility of the system at order 

In this appendix, we will show that equations given by the identification of the coeffi- 
cient (x — Sj)° given before do not add any new constraints. Indeed, from the equation 
(EH]) we get: 

v N ) 

.^rr'f \ Id i V" (^ + rw^)K(x ,s j ) 
+rw i K (x ,Si) = — h J — (C.l) 

Xq &i .,■ &i Sj 

A possible worry is that one can eliminate the term K'(xq, si) in the previous equation 
by multiplying on the left by rwj. Then one can obtain a closed system of equations 



5N 



for the K(xq, Sj)'s which is different from the one we presented before. We will prove 
here that this does not happen. 

Let's first project the equation by multiplying on the left by fwi . The first term 
disappear because of the relation (IB. 171) . Therefore we find: 

1™/a-<(, , Si) = -5*1 + 1 v W + ftg/)^.^ (c . 2) 

iV x - Si iV ^-7" - s^- 

From there we can extract rw* K' {x Gl Sj) and put it back into the first equation. We 
get a closed relation for the K(xq, Sj): 



U7i* T \ - (ni?/ + fw j t )K(xo, Sj) 



?■„ — « ■ AT 



x - Si N ^ Si - sj 

3Ti 



T T ^ ^ + fw/)K(x , Sj ) 



N(x -s^ N Si - sj 



3& 

Some simplifications occur and we find: 



T - - t m IT -->t 

jj-rw. 



T V- (-R-rWi )K(x 0l Sj) 



x - Si N <f-f S* - s, 

We can project this relation into components. Note that the first d<i — 2 lines do 
not give any new information since they are the same as (1B.14I) : 

(-v{{si) + —) K ^ s ^-^ J2 K{xo,Sj)p,q +K ^ s *w>* = 

\ N Si — Sj N Si — Sj Xq — Si 

(C.5) 

The projection on the last line is interesting: 



r 

N 



T . 1 \ 

Vl{si) + — — K ( x o, Si)d 2 -i, q 
iy j^i Si S 'J J 

d2-l / u \ 

Nt d2 k=Q \ j7 ,i S i~ S 3 J 

N d2 -l,q id2 f T 2 K(X , Sj) d 2 -l,g 

a. _ 



'/ 



+ AT2 



;Eo — Sj A?" 2 ^—f. Si — s 



/t-1 <* 2 -l / 



a2 fc=0 
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(C.6) 



In the last term, the case k = d 2 — 1 cancels the previous one: 



s i)d 2 -l,. 



rp <*2-l / 

Nt r ^ ' ^ 

r°d,2—l,q 1 



d 2 k=0 \j# Si S 3 



h + S k ,0Si \ K(x ,Si) kt q 
d 2 ~2 



n' 



q . i 



T 



x - Si 



EE 



(C.7) 



Nt d2 — 8i- 8j 



K(x , s 



j)k,q 



We can put this system into a matrix d 2 N x d 2 N form. We introduce: 

/ \ 

/ K(x ,s 1 ) hq \ 

K(x , s N ) 1>q 
K{x ,s 1 ) 2 , q 

K(x , s N ) 2>q 



K[X(j, Si)rf 2 _i i9 

\K(x , S N ) d2 -i, q ) 



9q 



1 



X -Sl 



x -s N 



tdr, {XO—Sl) 



—Uq,N 

\td 2 (^o-sjv)/ 



9(h-i = 



and 



f-B Id 
-B Id 



V 



\ 





: Id : 

... -B Id 

\ Co Ci ... w 2 _2 c<i 2 _i y 

where all the matrices C&, 7<i and are of size N x N given by : 



(C. 



(C.9) 



Ci. 



(Cfc)i ; i 
(C/c)i, 



iVt, 



Vtt- - t k + <5 fc , Si 



T Mfc.i 



(CIO) 
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and a special case for Cd 2 - 



i- 



r Yl 1 



Cd 2 -\ = j-; diag(-V[(s 



N 



Si S j 



b d 2 -l 



(C.11) 



Then the two sets of equations are equivalent to the matrix equation: 



Wq : VK q = g q 



(C.12) 



The main problem now is that this system is different from the one we got before 
and thus the system seems overdetermined. Fortunately, the matrices, Ai and V and 
the vectors f q and g q have the same d,2 — 2 first lines. But the main difference is that 
f q has double poles in x = s, on the last line whereas g q has simple poles. Writing: 



VM~ l f q 



(C.13) 



And looking at the last component gives that the two matrix systems are compatible 
if and only if the last column of is null. From the form (companion-like) of the 

matrices we have: It is easy computation to see that the last column can be split in 
blocks of size m : 

fx . 

Mr 1 = : 



X M \ 



(C.14) 



where: 



B k (B + B 1 B + --- + B^xB^- 1 ) 1 



(C.15) 



so that: Multiplying on the left by V gives that the two systems are compatible if and 
only if: 



C + C X B + 



Ofa-iS*- 1 = o 



(C.16) 



Note that we have also the determinantal identity: 



det(M) = det 



(-X 




\^4o 



Id 
-X Id 



Ax 



0\ 




Id ': 
-X Id 
A d -i A d 



-l) Nd det I ^AiX* 



(C.17) 



,i=0 



valid for every matrices Ai and every matrix X of size N x N. 

This condition seems rather completely non trivial from the definition of B and of 
Ck- But we will show in the following that from the definition of the Uk/s, the last 
equation is automatically satisfied. Note then that because of (IC.17P we obtain at the 
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same time that V is not invertible and thus only (1B.26|) gives us the K(xq, Sj)'s. Let's 
now prove that we have always: 



C + C X B + ■ ■ ■ + C da -\B d *- x = 
From the definition of the j's we have: 



(C.18) 



M d2 ~ 1 N d2 ~ 2 ( T 1 \ 

fc=0 fc=0 V r^i 4 r / 

, „ _/. *i *r 



We obtain then that: 



'd 2 -2 



T 



. fc=0 



'J 



(C.20) 



ivfcd 2 fc=o \r^i 



u Kr {B k \j + u Ki {B k 



Sj, Sy 



t k (B k ) i:j + sA, 



Note then that: 



N 



S a S if A- 



so that: 

^2-2 



r 



rf 2 -2 



. fc=0 



k=0 V 



NV{( Si ) 
T 



NV{( Sl ) 
T 



T 



Si s r 



~(B k )r,j — tk(B k )ij + SiS, 



■^u d2 _ 2 ,i(B° 

td 2 



(C.21) 



Then we need to add the last term 



N 



From the definition of Udo-i 2 we have also: 



N 



T ~ 



Ud 2 -2,i = ^~](B t )i,kUd 2 -l,k + Jyt 



k=l 



d 2 -l 



so that: 



1 T . 2T 2 ^ 1 Tt d ^ 
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and thus: 

(C.22) 

proving that only (IB. 26ft gives the correct formula to determine all the K (xq, Sj)'s. 



D Example: Only one root s 

In order to illustrate our algorithm, we present here the case when there is only one 
root s (we omit the subscript S\). In this case, almost all the computation can be 
carried out analytically and in particular we show how to obtain W^\x) and Fx in 
this context. The results presented here are direct applications of the formula presented 
in the article. First we have from (16.141) : 

dz-k-l 



B = V{{s) 



T 



Then we have: 



d 2 p-1 



p=0 



N(x 



p=0 q=0 



The action S is given by: 

S(s,s,A,u) = Vx(s) + V 2 {s) - ss 



The extremum is obtained for u = 1, A — 1, V{(s) = s and V^(s) = s 
The Hessian matrix is given by: 

fVf(s) 
-1 



"H(s, s, A, u) 
At the extremum it gives: 



-1 

V"(s) 












n 



|extr 



-i 



V o 



-l 






o \ 








T 


T 


NA 2 


N 


T 
N 


o / 





\ 



VZ{V{{s)) 












T 



_ T 

N, N 



which gives: 



/ 



Vi'(s)V{'{s)-l V-(s)V{'(s)-l 

1 V(' (8) 

V£(5) V{'(s)-1 V 2 " (a) V/' (s) - 1 










\ 










N 
T 


JV 
7M 2 / 



(D.l) 



(D.3) 
V£(V/( a )). 

(D.4) 



(D.5) 



(D.6) 
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so that at the extremum we get: 

/ Vj'(V((s)) 
I V4'(V,'M)V,"(s)-l 

H 1 |cxtr = 



2 (V{(s))V{'(s)-l V'i(V{(s))V"(s)-l 

V 2 "(V{(s))V{'(s)-l V>>(V{(s))V{'(s)-l 







\ 




-f 

N N 



(D.7) 



Hence we have with f!6.45p : 

ds 



dVi(x) (V 2 "(V{(s))V{'(s) - 1) (x - s) 2 
giving the formula for W 2 (%i > #2) : 



(D.8) 



W 2 {0 \x u x 2 ) 



Vi\V{{s)) 



(D.9) 



N {V 2 "{V>{s))V{>{s) - 1) ( Xl - sY(x 2 - s) 2 

We can now compute the kernels G(xq, x) and K(xq, Sj). The matrix M. (size d 2 x d 2 
is given by: 



/ 



M 



-V{(s) 






1 

-V{(s) 





1 



-V{(s) 







1 



\ 



(D.10) 



and the linear equation giving the i^(xo, s) reads for q = d 2 — 1: 

( K(x ,s) , d2 -x \ / \ 

^(^0, s )l,rf 2 -l 



M 



K(x Q , s) d2 - 242 -i 




1 



(D.ll) 



\(*o-s) 2 / 

The expression of A4.~ l is quite complicated. But fortunately, the last column 
(which is the only one we need for the computation of W 7 ^) is the most simple: 

v{{ s y 



(AT 1 ) 



-v('(s) e 

Vfe=o 



(D.12) 



T 

A? 



p=0 



But since it^i — ■»? £ ^d2+p+i(W( s )) P ' we hud easily that: 



NV{(s) p 



T{V»{s)V»{V{{s))-l) 



(D.13) 
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We can compute Lf[°\x, y\ x) by derivation of U^ v> (x, y) relatively to Vi(x) with the 



(0), 



formula 



U[°'(x,y;x) 



dV 1 (z) 



u[°\x,y; z). We find: 



T 



\q>(v>(s)) 



&i p-i 



N(V>i(V{(s))V{>(s)-l)(x-sr 

d 2 p-1 



T 



N(x - s) 



p=0 g=0 

dV{(s) 



T 



p=0 q=0 



^i(x) 



TV - 1) (x - a) 



p=0 q=0 



d 2 p-1 



EE^^^r 1 



N (* ~ s ) P=0 9 =0 

V{\s)V^{V{{s)) 



(x — s) 2 



- 1) (x - *) 



(D.14) 

If we project along the k th power of y, it gives: 



£/$(*; s) 



r 



r— 1 1 



(D.15) 



Therefore we have with (I7.16P : 
/ U$(x) \ 



Res K\x, z) 



( U[% ]Z ) \ ( U®{z) \ 



T 



(D.16) 



We want to compute W^\x) = j-U^^^x) that is to say that we only need the 



i rrW 



t,i. 



last column of the previous vector. If we write it in components we get: 



da-l 



fc=0 



t d M 1] {x) = Rff E K(x, z) Kd2 ^ ( U§(z; z) - *-d z U^{z) 



N 



r(0). 



(D.17) 



The contribution from Uq°1(z) is easy to compute since it only involves a double 
pole at x — s. Regarding the fact that K'(xq, s) P)q are null for p < d 2 — 1, there is only 
one term involved. An easy computation leads to: 



N 



d2 ~ 1 N~ 1 
Res E K (x,z)k,d2-idzU^°l(z) = —t d2 

a-m — ' i X — S 



k=0 



(D.18) 
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The contribution from U^.{z\z) is technically more involved. Indeed, since 
U{l(z\ z) has poles up to order 4, we need to compute the derivatives up to K"'(x , s). 
In particular, we only need the last column of the matrices. Since m = 1, the formulas 
of the derivatives simplify a lot: 



/ U ,i \ 



w 



td 2 



V = ^~ 

td 2 



( -k + s\ 

-tl 
\ -td 2 -l ) 



V 1 



d-2 





V o y 



V 2 = 



(D.19) 



which gives: 



K(x ,s) 

P,d2—l 

K'(x ,s) d2 _ 1>q 
K"(x , s) P)d2 _i 



-t 



1 



d2 (v 1 "(3)y a "(v 1 '(s))-i) (^o-s) 2 

N 2 u q>1 



TH d2 (x -s) 

^ (Sp td2 - 2 K'(xo, s) d2 - hd2 -i - V{'(s)K(x , s) pA2 -i) 
2NHp ^- 2 + r. (P<d 2 - 1) 



T 2 {x -s) ^ T(V 1 "(s)y 2 "(V 1 '(s))-l)(x -s) 2 



(D.20) 



Then ( 1B.29|) reduces for p < d,2 — 1 to: 



T_ is in 

3N 



K"'(x ,s) P!d2 . 



■\V[(s)K"(x , s) Pyd2 „i + \K"(xq, s) p+M2 _i 
V{"{») 



K(x , s) Pt d 2 -i 



(D.21) 



The last component is a bit more sophisticated. We can obtain it by multiplying 
(IB.29j) by rw* in order to get rid of K"(xq, s): 



6iV 



fWK"\x Q ,s) 



V"{s)tv? + — fV 1 K'(x Q , s 



N 

fw t K{x Q , s) Pid2 _i 



rw 



(Xq - S) 



(D.22) 



which gives in components: 



-^ 2 K"'(x , s) rf2 _i irf2 _i 
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k=0 



td 2 



(D.23) 

Now we can compute Wi'\xq). First we can compute the term coming from ^rrjp of 
f lLTTo]) in flDTfj) : 



(1) 



T 



cfo — 1 d 2 — k— 1 



Ar(^(y/( S ))^( S )-1) 



_i ^"(xq, s) k . 



d 2 -l 



k=0 r=0 
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j cfe-1 da-k-X 



M d2 t d2 vi'(s)\q»(v{(s)) 



~ T (V»(V{(s))V{>(s) - 1) (so - s) 2 (^(^'(^(s)) - l) 2 (x - s) 2 

(D.24) 

Then we observe that we have to compute: 

d'2 — l d2—k—X 



TV»{V'{s)) 

)-l) 

(E E W fc+ i(^(s)) r A^(x , S ) M2 _i 



(2) = ^(^^(^.^ E E ww^v. 



6iV(y/'( s )y/(T//( s )) - 1) 

+t d2 K'"(xo, s)d 2 -i,d 2 -xj 

TVf(V{(s)) 
W(V{>(s)V»(V{(s)) - 1) 



d 2 -2 d 2 -fc-l 



fc=0 r=0 

d 2 -2 



(E E w fc+ i(^( S )) r A^(x , S ) M2 _i 



-j^K (x , s) Ma _i — ^ (s)A (x , s) d2 _i id2 _i 

fc=0 



E u k,iK(x , s) k42 -i - 



s) 3 



k=0 

Vx( s ) K '( x o,s) d2 -x,d 2 -x 



TVZ{V{{s)) ( 6Nt d2 



6N(V{'(s)\q'(V{(s)) - 1) V T 



s) 3 



t d2 1/ 1 // (s)A"(x ,s) d2 _i id2 _ 



i 



(D.25) 



iV d2 ~ 1 £ 



where we have used the definition of w^i to cancel the terms in K"'(xq, s). Eventually 
we observe that: 

d 2 -x 

^ Mfc,iA'(xo, s)m2-i 

k=0 

rp7 d 2 -X 
1 l d 2 



E WW)' 



JV(VY^(V?( a ))-l)(so-a) a A 

<fe — 1 d 2 — 1 

m , m ,, (ww)-i)(xo- S ) 2 £ S Wfc+1 (v?(a))r+ * 
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(D.26) 



Ttd 2 


N(V{' 


[s)V»(V{(s))-l)(x - 






Tta 2 




N(V{' 


[s)V»{V>{s))-l){x Q - 


s) 2 




Ti d2 Vi\V{{s)) 





N(vf(s)v»(v;(s))-i)(x -sy 



d 2 -id 2 -i 

k=0 p=k 
p=0 k=0 



(D.27) 



Therefore we find: 



(2) 



Nt d2 Vf(s)V»(Vl(s)) 



+ 



t d2 vf(s) (\q'(v{(s))y 



T(V('( S )\q'(V!(s))-l)(x - S ) ' 2(V{'(s)V^V{(s))-l) 2 (x -sy 

wrcw) (D . 28) 



(V{>(s)\q>(V{(s))-l)(x -sr 

And eventually putting (1) and (2) together and the contribution from d z U^l{z), we 
get to: 



W[ l \x ) 



N 



T(V^V{(s))V('(s)-l)(x -s) 
V['{s)V-{V{{s)) 



2(v?( s )v»(v{(s))-iy(x - s y 



vf(s) W(v{(s))y 



NV»(s)V»(y{{s)) 



2 (V;>(s)Vl(V{(s)) - l) 2 (x - s) 2 T(V;i(s)V'i(V{(s)) - l)(x - s) 
V^V[{s)) , N 1 

t d 2 



(D.29) 



(V{>(s)V>>(V{(s))-l)(x -s)z T a2 x -s 

(v;'(s)v^v;(s))-i)(x - s f 

V{'{s)Vj"{V{{s)) + V{"{s) {V^V{{s))f 
2{V(\s)Vi'{V{{s))-lf{x - S y 



i.e.: 



W^ixo) 



V{' W (V{ (s))+V{" (s) ( VJ_ (V{(s))) ' 
2{V{'{s)V 2 "{V{(s))-iy{x -s)* 



(D.30) 

Then a straightforward computation gives F\ defined by a ^ x ^ F\ = (x) with: 



(V{'(s)V 2 "(V{(s))-l)(x -s)3 



F 1 = -- In (V?(s)V?(V{(s)) - 1) = -\ In (V^V^s) - 1) 



(D.31) 



where F Q is given by —S which reduces when N = 1 to: 



F = L (- Vl (s) - V 2 {V{{s)) + sV{{s)) = ^ (-Vt(s) - V 2 (s) + s~s) 



(D.32) 
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and which satisfies q^F = w[ Q \x) = N ^_ s) ■ 

E Computation of W^\x) in the general case 

We have: 

d 2 



[^(B)] i = E**( B( * > )< = a < 



fc=0 



where 



The relation 



implies 



where: 



A' 



3=1 

SV x [x) 

Y. Zi,jUj,k = V i>k 



dV{( Sj ) 



d(jV>(B)] l -s l 

(ISa 



The equation to solve is: 

UW(x,y) 



J2 — -V((x) + y 



T 

N ^— ' x — Si 



i -j 



vm-n^) i + 



+U^\x,x,y) = -P^\x,y) 



1,3 
(1), 



wP(x 



where 

r(0) 



SU^(x,y) _ 5 



U\ '(x,z,y) 



SVi(z) SV^z) 



N Z^ l >- 



vj{y) - v^jj) i 
y Sj % Si 



x-Vi 



Since we must perform an expansion in y, we write: 

d 2 -l 



Y, E j,k( s j^n)y k 

k=0 



and 



y~ s 3 



k=0 
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with 



$i,k = T7 E A i: jE jik (Sj, t n ) 



Since: 



EjtP—li.Sjytn) SjEjp(^Sj,t n ) 



td 2 

t p VI < p = d 2 — 1 
to - V^Sj) 



we obtain: 



^4^2-1 _ td 2 

%,p-i - E B ^k, P = t p Vl<p = d 2 -l 



E Bi^k,o 



t - [Vi{B)]i = t - Si 



We have the following properties: 



■p 



d 2 —p—1 

E < 

k=0 



k+p+1 



B 



(k) 



We calculate: 



E B«i™^ 



p=0 <j,i 



Consequently, 



E tk+p+iB 



(k) 



d.2—1 ,„ d 2 -p— 1 

E E^i,^ 

p=0 q,i ^ fc=0 

EEw^ E U-" 1 ' 

p=0 g,i u=p+l 
d2 it— 1 



EE*- E^ 

A~ _ d[VJ(B)] t 

/ 



q,i g(u-p-l) 



u=0 



dSj 
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We note that: 



P =o q>i dsk ds i 
d ^ d [ B P] d$ ijP 

p=0 q,i 



p=l q,i 



dsk dsj 



ds k 



1 dsj 



dsk ds^ 



2^ 2^ ds k ds j ^ ^ 

p=0 q,i J p=0 q,t 
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r ' q dSk dSj 



(E.16) 



so that 



dsjdsk 



d 2 -l 



E 

+ EE> B ' 
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ds 7 Sfc 



g?s 7 - dsk 
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Similarly, 



*2-l 



p=0 q,i u J ' p=0 



d[Vj{B)l 
dV{( Sj ) 



Vr,j 
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Now we have 

d 2 -l 

Ui°\x,z,y)=J2y P Yl 

p=0 i,k 

where 



A k $ iiP —^— — ^ + $i,p^i, fc7 — - 



^1(2 



X - Si (2 - s fc ) 5 



X - 



- s fc ) s 



Then, u[°\x, x, y) is obtained by changing z — > x and y — > y = — jj4: with the 
convention that y is written at the left of all x dependent terms. 

1 



(E.19) 
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1 1 1 
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Also we write 



d 2 -l 
p=0 



x) 



(E.22) 



so that 



d 2 -2 



UW(x,y)y = y d2 U d2 ^{x) + ^y f 

p=0 



U p (x) + dU ^ {x) 
ax 



+ 



dU x (x) 
dx 
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The system to be solved is: 



d 2 -2 



y^U^x) + J2 



p=0 



U p (x) + 



dU p+1 (x) 



dx 



d 2 -l 
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Clearly 



U d2 . l (x)=t d2 Wi 1 \x) , jj-^u- 



T d 
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and 



We have 



Wi 1] (x) 
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X-Si (x-s k ) 2 
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1 1 
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td 2 
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x - Si(x - s k )' 



(x - s^ 2 (x - s k y 



The term involving the derivatives is easy to compute. Indeed, it gives: 



N 



2^ K '( X ^ S i)k,d 2 -l U k; 



(E.30) 



2 i k=0 



TL 



Now only the term with k = d,2 — 1 contributes from flB.30j) . Moreover, Ud 2 -i,i - 
and from ( IB. 321) we have: 



K'(xq, Si)d 2 -i,d 2 -i 



N 



T(x - Si 



K(x , s 



j)d 2 -l,d 2 -l 



Si Sj 



so that the term involving the derivative is: 



(E.31) 
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N 

N 1 
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-2E 
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Note in particular that it contributes with a simple pole at x = Sj with residue 



where the kernel K(xq, x) satisfies 



V 1 

Ndx 



Id 



K(x 0l x) = G(x , x) = 1 

z — ' X — Si Xq — X 
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the matrices ai and Id are d,2 x d,2 matrices. We expand around s,-: 



1 1 
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X Si (x Sfc) {x Sj / 
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(E.35) 



Consequently, from the expansion of K(xq,x) around x = Sj\ 



<h 



with 



j p=i 



(E.36) 



(E.37) 



Si Sj 
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The expansion of X>* around x = Sj is: 
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where we have noted: 
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We also expand the kernels: 



G(x , x) 



EE 



- stAP 



[x — s 



+ E 



-Id (E.44) 
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We define the matrices: 
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\$ ii0 'I',. i ••• $i,d 2 -2 'I',../, :/ 
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We note that 
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MjCti 
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T 

N 



T 

N 



The c?2 — 1 first equations tell that: 
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We proved before that this is equivalent to 

dW{B)]t 



k 



r,(l,d 2 )" 



(x - Sj) 2 



(E.50) 



so that: 



and 
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E 
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The other equations tell 
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Finally, 
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We need only to look at the last equation: 

■y d 2 d 2 . d 2 

da 



p=l 



p=l i 



p=l i 



p=l 



(ar - Sj) 3 



(E.55) 



We now show that in Ud 2 -i, if we replace the terms in k'"^ d ^ and k"^ p d ^ by their value, 
then the terms in k\ , d s disappear so that only the terms in ki^ p ^ 2 ) which are known 
remain. 

First Step 



9 ^ y? jJ»-i' e j,(p,«fe) s -s- 



i 
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where, for p = 1, 
Step 2 



[A^ ijP _ 2 ] p _ x Aj[y 2 '(S)]j = A jSj = H 
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d2 
i p=l 

(E.59) 
Step 3 

From the indices properties of Ch k we obtain: 



d 2 
p=l 



.)■:> 



T 



(x - Si y " a2 N(x - Sj ) 2 

~ \| \H-j,jki,(p,d2) + 2"Hi,jfcj,(p,d 2 )) 



'■j,(p,d2) 



d 2 



— s,, 



(Si - Sj) 2 



rf 2 

i,j p=l 



i,p-l T *j,p-l 
Sj 



fcj,(p,d 2 ) 



(E.60) 



Consequently, 




(E.61) 



Note in particular that the simple pole in ( 1E.60I) is canceled by the one coming from 
flK32|) 



We now calculate Bj, we have 
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- E j^7f z J* h B '-\* ( - % E *t*.«w 
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The bracket is 
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Consequently, 
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We have: 
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We note that 
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Consequently, we find: 
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We have not been able to sum over p the last term. This result is the same found from 
the differential operator y = — 4- method. We have 
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and write: 
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